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ABSTRACT 

We study the dependence of angular two-point correlation functions on stellar mass (M*) and 
specific star formation rate (sSFR) of M* > 10 1 o Mq galaxies at z ~ 1. The data from UKIDSS 
DXS and CFHTLS covering 8.2 deg 2 sample scales larger than 100 /i _1 Mpc at z ~ 1, allowing us 
to investigate the correlation between clustering, M*, and star formation through halo modeling. 
Based on halo occupation distributions (HODs) of M* threshold samples, we derive HODs for 
M* binned galaxies, and then calculate the M*/Mh a io ratio. The ratio for central galaxies shows 
a peak at Mhaio ~ 10 12 /i _1 Mq, and satellites predominantly contribute to the total stellar mass 
in cluster environments with M*/Mh a i 0 values of 0.01-0.02. Using star-forming galaxies split by 
sSFR, we find that main sequence galaxies (logsSFR/yr _1 ~ —9) are mainly central galaxies in 
~ 10 12 5 /i _1 Mq haloes with the lowest clustering amplitude, while lower sSFR galaxies consist of 
a mixture of both central and satellite galaxies where those with the lowest M * are predominantly 
satellites influenced by their environment. Considering the lowest Mh a i 0 samples in each M* bin, 
massive central galaxies reside in more massive haloes with lower sSFRs than low mass ones, 
indicating star-forming central galaxies evolve from a low M*-high sSFR to a high M* low 
sSFR regime. We also find that the most rapidly star-forming galaxies (logsSFR/yr _1 > —8.5) 
are in more massive haloes than main sequence ones, possibly implying galaxy mergers in dense 
environments are driving the active star formation. These results support the conclusion that the 
majority of star-forming galaxies follow secular evolution through the sustained but decreasing 
formation of stars. 

Subject headings: galaxies: evolution — galaxies: halos — large scale structure of universe 


1. INTRODUCTION 

It is expected that small structures merge to 
form more massive ones in the lambda cold dark 
matter (ACDM) paradigm. Therefore small dark 
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matter haloes are the seeds for larger structures. 
Galaxies form in these dark matter haloes through 
the binding of baryo ns and the cooling of gas 


(White & Rees 1197811 . Since galaxies evolve in 


their host haloes, the distribution and evolution of 
galaxies are tightly relate d to their host dark mat¬ 


ter haloes (Baugh 1200 6). In the context of hier¬ 


archical structure formation models, massive dark 
matter haloes can contain many galaxies with a 
wide range in mass. Furthermore, the population 
of member galaxies depends on the properties of 
their host dark matter halo, because the potential 
well of the dark matter halo affects the properties 
of galaxies within it. 

A popular method for measuring the distribu- 
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tion of galaxies is the two-point correlation func¬ 
tion, which describes the excess probability of a 
galaxy pair over a ran dom distribution on specific 


minous galaxies arc found in more massive haloes 


scales ( Peebles 1980l) . Recently, wide and deep 


(Zhene: et al. 200?]: 

Ross & Brunner 2009a; Wake et al. 

2011; Zehavi et al. 

2011; CouDon et al. 2012). 


surveys have provided opportunities to study the 
dependence of the clustering of galaxies on their 


- - — — to" ~to~^ — --— , — 

ity, stellar mass and population (Norberg et al. 

2001 

, 20021 Zehavi et al. 2002, 

20051; Coil et al. 

2008 

: Ross & Brunner 2009al 

Loh et al. 201C& 

Ross et al. 2010: Zehavi et al. 

20111. As a re- 


suit, it is known that redder, brighter or more 
massive galaxies are more strongly clustered than 
those having opposite properties. 

For high redshift galaxies, one of the most 
efficient selection methods is to use their ob- 


However, the HOD framework fits the number 
density of galaxies and their clustering simulta¬ 
neously so we require a large number of galaxies 
to sufficiently constrain the parameters of both. 
Also, to avoid the effects of cosmic variance, these 
galaxies need to be mapped over a large area of 
sky. The lack of large, sensitive near-IR detec¬ 
tors has prevented identification of a large num¬ 
ber of galaxies at z > 1, where the bulk of stel¬ 
lar emission is observed in a near-IR regime. So 
far, previous work on stellar mass limited galaxies 
at high redshifts has been based on survey data 
with areas from a few hundreds arcmin 2 to <1.7 


2000: Im et al. 


20041; Kong et §1. 


McCracken et al. 


Merson et al. 


served color. Thus, many previous studies have 

deg 2 (Meneux et al. 

200a 

2009; 

Foucaud et al. 

applied various color cuts to select high red- 

201C: Hartley et al. 

2010; 

Furusawa et al. 2010; 

shift galaxies such as Extremely Red Obiects 

Wake et al. 2011; Hartley et al. 

20130. Further- 


more, this small surveyed area makes it difficult 
to measure reliable clustering strengths on larger 


2013; Kim et al. 112014), BzK galaxies (DacrLetal 


_ _ _ Daddi et al _ _ 

2002: Roche et al. I 2002|,~ 2003; Yanet_al^ [2004; 

Brown et al. 12005: Kon g et al. 1120061: Gonzalez-Perez et akcldes (a few tens /i -1 Mpc), where the distribution 
2009; Kong^toil. 2009^ Kim et al. 20 lF Palamara et al. b f dark matter haloes is imprinted, since these sur¬ 
veys have covered at most ~ 50 /i -1 Mpc at z ~ 1 
on a side. Therefore it is important to perform this 
analysis with homogeneous galaxy samples drawn 
from a wide-area near-IR survey. 

In terms of galaxy evolution, differ ent mecha¬ 
nisms play key roles at different epochs. IPeng et ~ 


2006; 

201oF 


ies (Franx et al 


2013)_ and 
2003 


lHartlev et al. 

Hanami et al. 

Distant Red 


Grazian et al. 


Galax- 

20001 : 


Fouca ud et al. 200'mQuadri et al. 12008tlGuo fc White 


20091) . Their clustering properties also show sim¬ 
ilar trends to those of low redshift galaxies. Al¬ 
though color selection is efficient in isolating galax¬ 
ies in a specific redshift range, a simple color cut 
often extracts a mix of galaxies with different 
properties and redshift. Hence to correctly trace 
galaxy clustering it is necessary to measure it from 
a sample containing galaxies with well defined in¬ 
trinsic properties and a narrower range in redshift. 

The halo occupation distribution (HOD) frame¬ 
work makes it possible to interpret the galaxy clus- 
tering in relation to their host dark matter haloes 
(Ijing. Mo fe Boerner 1998 : Benson et al. 2000: 


( 2010t l proposed that mass quenching is impor¬ 


tant for all galaxies, but environment quench¬ 
ing dominates at low redshift and at lower 
masses. Also, there have been many results 
about the stellar mass function of passive galax¬ 
ies from wide field optical-near-IR datasets, 
showing relatively mild evolution of the most 


Ma fe Fry I 20001: Peacock fe Smith 2000t Seliak 

2001)1: IScoccimarro et al 112001 : Berlind fc Weinberg 


low mass ones (Drorv et al. 2009c Ilbert et al. 

201C: Bezanson et al. 2012; Ilbert et al. 

2013; 

Moustakas et al. 2013; Muzzin et al. 2013 

Tinker et al. 

2013: Tomczak et al. 20141), which is consistent 

with the model in Peng et al. 

201 Oil. In addi- 


Goorav fe Shethl l2002h . The HOD quan- 


tion, the evolution of the luminosity function of 
Luminous Red Galaxies (LRGs) foll ows a pass ive 
evolution model at z < 0.6 ( Wake et al. 200fith 


tifies the mean number of central or satellite 


galax i es in a given halo mass (IKravtsov et al. 
2004 : [Zheng et al. l2005h . Based on the halo 


model analysis with multiwavelength datasets, 
many authors have reported that massive or lu- 


IScoville et al. I ( 2013l l reported that the fraction of 
early type galaxies increases from 30% at 2 ~ 1.1 
to 80% at z ~ 0.2 in the densest regions, but from 
30% to only 50% in low density regions. These re¬ 
sults suggest that the dependence of galaxy prop- 
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erties on their host haloes at z ~ 1 is different 
from that in the local universe. 

Furthermore, the relation between galaxy prop¬ 
erties and their environment a t z ~ 1 is still con ¬ 
troversial. Elbaz et al. ( 2007 ) and lCooper et ah 


(120081 ) found a reversed relation between star for¬ 
mation rate (SFR) and environment at z ~ 1, 
meaning a higher SF R was observed in th e high¬ 


est density regions. IScoville et al. ( 20131 ) found 


the evolution of the relation as a function of red- 
shift, and recover a weak or no dependence of 
SFR on environment at z ~ 1. On the other 


In § [2l we briefly describe each survey and how 
the catalog was generated. In § QJ| we note the 
methods applied to select samples, to measure 
clustering and to model the HOD. The depen¬ 
dence of galaxy clustering on stellar mass and the 
stellar mass to halo mass relation are described 
in § |4] We present the dependence on sSFRs in 
§ EJ We also relate these results to the evolution 
of galaxies in § El and finally summarize this work 
in § 0 Throughout this paper, M* indicates a 
stellar mass of galaxies and M^aio means a dark 
matter halo mass. We assume a flat ACDM cos- 


hand. CooDer et al 

(2010l), Chuter et al. 

20 111) 

mology: Cl m = 0.27, ug = 0.8, Ho = lOO/i km s 1 

and Quadri et al. 

(2012h reported that the color 

Mpc -1 with h = 0.71. The photometry is quoted 


(or S FR)-density relation persists out to z ~ 
1.5. iTinker et al. ( 20131) argued that the cen¬ 
tral galaxies in low mass haloes are likely to be 
star-forming galaxies at z ~ 1, and their evolution 
contributes to the observed change in the red se¬ 
quence. However, there are not many results con¬ 
necting galaxy properties with their host haloes 
at z ~ 1, especially detailed, statistical studies 
based on a l arge samp le from wide area surveys 
(>2-3deg 2 ). iMostek et al. I ( 2013 ) measured the 
clustering strength of galaxies with various criteria 
such as stellar mass, SFR and specific star forma¬ 
tion rate (sSFR) at z ~ 1 from a small area spec¬ 
troscopic survey. We re-address their work with a 
much larger photometric dataset making it possi¬ 
ble to split galaxies into finer sub-samples. 

In this work, we use wide and deep multiwave¬ 
length datasets with ugrizJ- and A'-bands based 
on UK Infrared Telescope (UKIRT) Infrared Deep 
Sky Survey (UKIDSS) Deep eXtragalactic Sur¬ 
vey (DXS) and Canada-France-Hawaii Telescope 
(CFHT) Legacy Survey (CFHTLS)-Wide. The 
catalog covers 8.2 deg 2 with the limit magnitude of 
Jab = 23.2. This is one of the best datasets to in¬ 
vestigate the clustering properties of homogeneous 
galaxy sub-samples and to minimize the influence 
of cosmic variance, thanks to its unique combina¬ 
tion of depth and area. Furthermore, these data 
allow us to link galaxies with various criteria to 
their host haloes separately. Using this catalog, we 
measure the angular two-point correlation func¬ 
tion of z ~ 1 galaxies split into several sub-samples 
based on stellar mass and sSFR. Additionally, we 
fit a halo model and measure the bias factor with 
the measured correlation function in order to link 
galaxies with their host dark matter haloes. 


in the AB system. 


2. DATA 


2 . 1 . UKIDSS DXS 


The DXS (A. C. Edge et al. 2015, in prepara¬ 
tion) is a sub-survey of the UKID SS w hich was 


perfo rmed from 2005 to 2012 (ILawrence et al. 


l2007t ) using the UKIRT. The DXS images were ob 
tained using the Wide Field Camera (WFCAM, 
Casali et al. 2007) composed of four Rockwell 
Hawaii-II 2Kx2K array detectors covering four 
13.7x13.7 arcmin 2 regions. Since WFCAM has 
a relatively large pixel scale as 0".4/pixel, a mi¬ 
crostepping technique has been applied so that a 
science image has O".2/pixel and avoids an under¬ 
sampled point spread function. 

The DXS maps ~35 deg 2 composed of 4 differ¬ 
ent 8.75 deg 2 patches (XMM-LSS, Elain-Nl, Lock- 
man Hole and SA22) with aimed depths of Jab = 
23.2 and Kab = 22.7 at a 5a point-source sensitiv¬ 
ity. The actual data show a 90% poi nt-so urc e com - 


pleteness at these magnitudes (IKim et ah 2011). 


The scientific goals of the survey are to determine 
the abundance of galaxy clusters at z > 1, to un¬ 
derstand the clustering of galaxies, and to inves¬ 
tigate the census of the luminosity density in star 
formation. 

In this study, we deal with the SA22 field cen¬ 
tered on a =22 h l7 m m s and S = +00°20'00"(J2000). 
In the whole surveyed area (3.4x2.6 deg 2 ), we per¬ 
form our analysis with images from the UKIDSS 
data release 9 (DR9) covering ~ 8 deg 2 , and one 
remaining WFCAM field (~0.7 deg 2 ) which was 
not released in DR9 comes from DR10. The area 
coverage corresponds to roughly 140 h -1 Mpc x 
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107 /i -1 Mpc at z = 1. The average seeing is ~ 0"8 
in both J and K. The photometric and astromet¬ 
ric solutions are based on the output from the 
standard pipeli ne, and are accurate to better than 
2% and ~0".05 ( Dve et al~ll2006 ; Lawrence et al. 


20071b 


2.2. CFHTLS 

The CFHTLf){j] is a set of deep and wide opti¬ 
cal surveys performed using the MegaCam cam¬ 
era mounted on the CFHT with ugriz filters. Of 
the three surveys that constitute the CFHTLS, 
we deal with the CFHTLS-Wide W4 field which 
covers the DXS SA22 area. The CFHTLS W4 
field covers 25 deg 1 2 with limiting magnitudes (50% 
completeness for point sources) of u ~ 26. 0, q ~ 
26.A r ~ 25.9, i ~ 25.7 and z ~ 24.6 (iGwvn I 
2012f) . For this work, we use CFHTLS images 


taken from the MegaPipe data pipeline at th e 
Canadian Astronomy Data Centre ( Gwvn I 2012T) . 
The image reduction procedure, as well as pho¬ 
tometric and astrometric calibrations are well de¬ 


scribed in Gwyn (12012 1. 


The UKIDSS DXS area is located in the south¬ 
east corner of the CFHTLS W4 held. Thus we 
extract only the sub-region that overlaps with the 
UKIDSS DXS held. We use the images in all the 
CFHT hlters, i.e., u,g,r,i, and z. 


2.3. Catalog 

Our main goals in this work are to measure the 
clustering of galaxies at z ~ 1 and to investigate 
how the clustering of galaxies correlates with the 
star formation activity and stellar mass of these 
galaxies. Therefore it is important to accurately 
determine colors of galaxies to perform spectral 
energy distribution (SED) hts to estimate galaxy 
properties as well as to determine their photomet¬ 
ric redshifts. In an attempt to improve the photo¬ 
metric accuracy, we generate a new catalog instead 
of using the released catalogs from the UKIDSS 
team via WFCAM Science Archive^. The new 
catalog of objects was constructed using the pro¬ 
cedure below. 

First, since huxes must be measured from the 
same region of galaxies at different bands to obtain 

1 http: //www.cfht.hawaii.edu/Science/CFHTLS/ 

2 http://surveys.roe.ac.uk/wsa/ 


accurate color, images are convolved through a 
Gaussian filtering to unify the FWHM. The worst 
seeing condition of our UKIRT dataset is used as 
a reference, which corresponds to a FWHM of l'/l. 
Then we resample all images to have the WF¬ 
CAM field of view with a four pointing observation 
and CFHT pixel sca le, using the Swarp software 


(Bertin et ah 200 2). Hence, we have 12 images 
for each band, each covering 0.8 deg 2 . We refer 
to each 0.8 deg 2 image as a tile. The coordinate 
offsets (AR.A. and Adecl.) between both surveys 
are 0"05 and 0'/04, respectively, which are accurate 
enough to align images. 


Second, we run the SExtractor software (Bertin & Arnouts 


19961 ) in dual mode on each tile of the covered 


area. The unconvolved J-band image is used for 
the detection and the measurement of AUTO mag¬ 
nitudes which we consider a proxy of total magni¬ 
tudes of an object. In order to measure the color of 
each object, 2" diameter apertures are used to de¬ 
rive aperture magnitudes on the Gaussian filtered 
images in all bands. 

Third, spurious objects such as cross-talk and 
objects on diffraction spikes are removed. It is well 
known that bright stars make cr oss-talk at specifi c 
positions on the WFCAM chip ( Dve et al. Il2006l ) . 
In addition, the unreliable regions of haloes or 
diffraction spikes of bright stars are a lso masked. 
We follow the algorithm in iKim et al. I ( 2011 ) for 
reducing spurious objects. The masking region 
for CFHTLS is taken from the Megapipe data 
pipeline. 

Finally, all sub-catalogs from different tiles are 
merged to create a single, master catalog. If ob¬ 
jects have been detected in an overlapping re¬ 
gion of different tiles, a 1" matching radius was 
applied to identify objects in common for both 
tiles and the weighted-mean of fluxes from differ¬ 
ent tiles are assigned as the final flux of uniquely 
identified objects. Additionally, the Galactic ex- 
tinction is corr e cted b ased on the dust map from 
Schlegel et al. I ( 19981) . In total, ~0.86 million ob¬ 
jects are found in the 8.2 deg 2 UKIDSS DXS area. 
Figure [l] shows the number counts of all galax¬ 
ies in the DXS area (solid line). Stars were ex¬ 
cluded using (g — z) and (z — AT) colors, which is 


similar to the BzK diagram (IDaddi et al. 1 12004; 


Oi et al. Il2014lh Figure [2] displays all detected ob¬ 
jects in a single tile (0.8 deg 2 ) with J < 22.5 as 
an illustration. The line indicates the criterion, 
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Fig. 1. - Number counts of galaxies in the DXS 
SA22 field. The solid line is for all gala xies in this 
work. Dashed line is the result from I. Tar vis et al.~l 
(12013 ). Our result is consistent with previous re¬ 
sults up to Jab ~ 22.5. 


(z — K) = 0.44(g — z) — 0.52, which we apply 
to distingui sh galaxies from st ars. The number 
counts from l.Tarvis et al. 1 ( 20H i are also indicated 
with a dashed fine in Figure [U Our result shows 
a relatively good agreement over the whole mag¬ 
nitude range. However, we may miss some galax¬ 
ies in the faint regime ( Jab > 22.5) due to the 
relatively shallow depth of CFHTLS and the as¬ 
sociated incompleteness. In fact, the distribution 
of objects shows a large scatter in their colors in 
the faintest magnitudes, where the magnitude of 
optical bands is close to the limiting magnitude, 
making it difficult to distinguish galaxies from 
the stellar locus. In addition, the 90% complete¬ 
ness for extended sources is ~0.3-0.5 magnitude 
brighter than that for point sources. We note that 
we will apply additional criteria of galaxy stellar 
mass (M* > 10 10 Mq) and photometric redshift 
(0.8 < z < 1.2) for our analysis, so that we select 
relatively bright galaxies which are well separated 
from the stellar locus so we expect negligible stel¬ 
lar contamination in our sample. Red dots in Fig¬ 
ure [2] are the objects satisfying our selection crite¬ 
ria. For the display purpose, we display a portion 
of full samples. 


Fig. 2.- The (z — K) versus (g — z) two-color 
diagram. We extract objects in a single tile (0.8 
deg 2 ) with J < 22.5 for a display purpose (black 
dots). The line indicates the criterion to distin¬ 
guish galaxies from stars. Stars fall below the solid 
line. Red dots are objects satisfying our selection 
criteria for the analysis in section 3.1. For the 
display purpose, ~10,000 objects are shown. 


3. METHODS 
3.1. Selection 

This section describes the basic measurement 
for galaxy properties and the criteria for our sam¬ 
ple selection with the multiwavelength catalog 
from it-band to A’-band corresponding to the rest- 
frame UV to near-IR for z ~ 1 galaxies. 


3.1.1. Photometric Redshift 


The redshift information is crucial to estimate 
galaxy properties and to investigate the clustering 
of galaxies in a specific redshift range. Using our 
multiwavelength catalog, we estimate the photo¬ 
metric redshift of each object. 


We used the Le Phare ( Arnouts et al. 19991: 


Ilbert et al. 20(1 fill software to derive photometric 
redshifts. We used 66 SED templates app l ied fo r 
the CFHTLS-Deep fields in Ilbert et al. I ( 2006h . 
The template s are based on E l l, Sbc , Scd and Irr 


spectra from 1 Coleman ct al. (19801) and a star- 


burst SED from iKinnev et al. I ( 19961) , and cover 
the wavelength range from rest-frame UV to near- 


5 































IR (see Ilbert et al. 2006 for more details). The 
Le Phare code produces offsets in magnitude in 
each band after running the code on a training set 
of galaxies for which spectroscopic redshifts are 
available. The application of the magnitude offset 
improves the photometric redshift accuracy, and 
as a training set, we used the VIMOS-VLT Deep 
Survey (VVDS) wide which is a spectroscopic sur¬ 
vey mapping 4 deg 2 of th e UKIDSS DXS SA22 
area for I < 22.5 o bjects ( Le Fevre et al. I [20051: 
Garilli et al. I l2008h . We select 3609 galaxies in 


the SA22 held having reliable spectroscopic red- 
shift information as a training set. First, we ran 
the Le Phare software for the cross-matched ob¬ 
jects between UKIDSS DXS and VVDS to calcu¬ 
late the magnitude offset for each band compared 
to those from templates, which are ~0.25 mag for 
J-band and less than 0.08 mag for the others. The 
large offset in J-ban d appear in other works (e.g., 


Ilbert et al. ( 20091 )). and can be understood as 


due to a template mismatch. In this case, the red- 
shift was fixed with that from VVDS. After find¬ 
ing and then applying these magnitude offsets to 
the objects in the catalog, the Le Phare software 
was run again for all detected objects. Figure [3] 
displays the comparison between measured photo¬ 
metric redshifts ( z p hot ) and VVDS spectroscopic 
redshifts ( z spec ) for the cross-matched galaxies. 
We find that the normalized median absolute de¬ 
viation of photometric redshift in Az/(1 + z) is 
~0.038. For galaxies with 0.8 < z spec < 1.2 that 
are studied in this work, the uncertainty is ~0.042. 
The fraction of outliers (|Az/(l + z)\ > 0.15) is 
< 5% for both cases. The dashed lines in Fig¬ 
ure [3] show where Az/(1 + z) = 0.15. Hereafter 
the term redshift (z) indicates the photometric re¬ 
sult ( z p hot )■ We use photometric redshifts only 
for the analysis. We also note that the magnitude 
offset is applied for the estimation of photomet¬ 
ric redshifts only. The magnitude offset improves 
the photometric redshift estimation by removing 
a systematic offset of Az =0.04. 

Finally, we apply a redshift cut of 0.8 < z < 1.2, 
and remove Galactic stars based on the gzK dia¬ 
gram as described in § 12.31 Note that no magni¬ 
tude cuts are applied when deriving photometric 
redshifts. In total 141,947 galaxies in this redshift 
range are used for estimating stellar masses and 
star formation activity. 

In order to check how efficiently the use of 



Fig. 3.- Photometric redshifts ( z p hot ) versus 
spectroscopic redshifts ( z spec ) for objects observed 
by the VVDS survey. The dashed lines indicate 
Az/(l + z) =0.15. 


photometric redshifts captures galaxies at 0.8 < 
< 1.2, we use a galaxy mock cata log from the 


GALFORM semi-analyt ical model ( Gole et al. 200rt 


Merson et al. 2013!. For this test, J-band mag¬ 
nitudes are randomly scattered in the photometric 
uncertainty ranges as a function of magnitude, and 
we selected galaxies satisfying J < 23.2, which is 
the magnitude limit in the observed catalog. We 
also randomly assign the redshift uncertainty to 
the mock galaxies. Then we compare the number 
of galaxies selected with modified redshifts to that 
with true redshifts. We find through this test that 
the use of photometric redshifts can recover 90% 
of galaxies with true redshifts at 0.8 < z < 1.2, 
while the interlopers (foreground or background 
galaxies) are about 13% among galaxies at 0.8 < 
z < 1.2. 


3.1.2. SED Fit 

In order to estimate stellar masses (M*) and 
SFRs of galaxies, we fit model templates of syn¬ 
thetic stellar populations to the multi -band pho; 


tomet r y fol lowing the algorithm of iLee et al. 


( 2010 . 2014 ). Here we briefly note the assump¬ 
tions made in this analysis. We u se SED templ ates 
from iBruzual fe Chariot I ( 20031) with a Chabrier 
( 20031 ) initial mass function. We assume a de- 


6 



































layed star formation history with an age (t) from 
0.1 Gyr to the age of Universe at the redshift 
of the galaxy in question and a star formation 
tim escale paramete r (t) f rom 0.1 to 10 Gyr. Also 
the Calzetti et al. ( 2000li dust attenuation curve 
is assumed for internal extinction. The reddening 
parameter of E(B — V) ranges from 0.0 to 1.5 
with a step size of 0.025. The metallicity was al¬ 
lowed to have values of 0.2, 0.4, 1.0 and 2.5 Z 0 . 
The SED fit returns the best fit parameters such 
as M *. SFR, age, r and E(B — V). The SFR is 
defined as the averaged one over recent 100 Myr, 
based on the reasoning of Lee et al. 1 (l2009blh Fig¬ 
ure |4] displays examples of SED fits for galaxies 
at 2 ~ 1. The solid line is the best fit SED, and 
points are observed fluxes at each band. The dot¬ 
ted line shows the SED in the top-left panel. For 
comparison purposes, we normalize the SEDs to 
the flux at observed frame 1.2 /iin of each object. 
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Fig. 4.— Examples of best fit SEDs (solid line) 
with observed fluxes (points) for galaxies at z ~ 1. 
Stellar mass and sSFR values of each galaxy are 
noted in each panel with units in M 0 and yr — , 
respectively. The dotted line is the passive galaxy 
SED in the top-left panel with normalizing to the 
flux at 1.2 /rm of each object. 

The left panel in Figure [5] displays sSFR versus 
M* of galaxies at 0.8 < z < 1.2 in the SA22 field. 
We find that the fit of a power-law (log sSFR = 
alog M* + /3) to the main sequence of star-forming 
galaxies has a slope of a = —0.33 ± 0.03 and an 
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intercept of /3 = —5.59±0.31 (red line). For the es¬ 
timation, we perform a double Gaussian fit to the 
sSFR distribution of galaxies in 5 different stellar 
mass bins from M* = 1O 1O M0 to M* = lO n ' 3 M 0 
with a width of 0.3 dex. Then the power-law is 
fit to the peak location of sSFRs for star-forming 
galaxies only. As a reference, the main sequence 
of star-forming galaxies at 0. 8 < 2 < 1.2 with the 

( 20073 ) 


power-law slope of -0.10 from lElbaz et al. 
is displayed with a blue line. In addition, we also 
plot the main sequence of star-forming galaxies 
at 0.5 < 2 < 1.5 taken from lWuvts et al.~T ( 201lh 
(green line), although the power-law slope of unity 
for the SFR-M* relation was assumed in their 
analys is. On the other hand, IWhitaker et al. I 
( 2012h suggested the significant evolution in the 
power-law slope of the SFR-M* relation with red- 
shifts. The relation at z = 1 derived by equa¬ 
tions in Whitaker et al. I (2012) is shown as an or¬ 
ange line with the slope of -0.43, which is similar 
to our estimate. We also note that our_relation 

d2012h . and 


is similar to that in IWhitaker et al. 
our sample includes mo re high mass galaxies than 
Whitaker et al. (2012) ■ This figure shows two loci 
of galaxies, one with logsSFR/yr -1 > —10, and 
another with logsSFR/yr -1 < —10. 

Based on Jab < 23.2 galaxies at 1.15 < 2 < 
1.2, the 80% percentile of the stellar mass distri¬ 
bution is M* ~ 10 9 ' 5 Mq. However, we may miss 
a fraction of the lower mass, passive galaxies due 
to the relatively shallow optical dataset. There¬ 
fore we use only galaxies more massive than M* = 
10 10 Mq for our analysis. The right panel in Fig¬ 
ure [5] shows the stellar mass of detected galaxies as 
a function of redshift with logsSFR/yr -1 < —10 
galaxies plotted in red and all other galaxies in 
black. The horizontal line is the stellar mass cut of 
M* = 1O 1O M 0 . Consequently, we extract 66,864 
galaxies in 0.8 < 2 < 1.2 with M* > 1O 1O M 0 for 
this work. 


3.2. Angular Two-point Correlation Function 

One of the simplest ways to measure the clus¬ 
tering of galaxies is the two-point correlation func¬ 
tion, which is the excess probability of finding a 
galaxy pair over a ra ndom distribution at a given 


scale ( Peebles 1980h . Here, we measure the an¬ 


gular two-point correlation function of galaxies at 
0.8 < 2 < 1.2, using the estimator introduced by 


7 

























































logM. (Mj) 


12: All 


logsSFR/yr'VlO 



0.8 0.9 1.0 1.1 1.2 

Aih-jl 


Fig. 5.— Left: the sSFR-stellar mass distribution of galaxies at 0.8 < 2 < 1.2 in the SA22 field. The red 
line indicates the sSFR of the main sequence galaxies in this work. Red points present the peak location 
for star-forming galaxies derived by a double Gaussian fit to galaxies in diff erent stell ar mass bins (see the 


text fo r m ore details). The blue, g reen and orange lines are results from lElbaz et al. I (|2007h , IWuvt.s et al. 


( 201 if ) and lWhitaker et al. ( 2012D . respectively. Right: the stellar masses of galaxies in 0.8 < z p hot < 1.2. 
Black and red dots are all and logsSFR/yr _1 < —10 objects, respectively. The cyan line is the stellar mass 
cut (M* = 10 10 Mq) applied in this work. For the display purpose, objects in 1 deg 2 of surveyed area are 
plotted. 


Landv fc Szalav 1 1 19931 ) : 


^obs(^) — 


DD(0) - 2 DR{9) + RR{6) 
RR(0) : 


(1) 


where DD is the number of galaxy pairs in 6 ± A 9 
in the observed data. In this work, the bin 
width is chosen as A log 0/degree = 0.15. DR 
and RR are the number of galaxy-random and 
random-random pairs, respectively. We gener¬ 
ated 30 different random point catalogs having the 
same areal coverage and angular mask as observed 
with each random point catalog containing a sim¬ 
ilar number of random points to that of observed 
galaxies. The errors on the two-point correlation 
functions and the covariance matrices were esti¬ 
mated by the Jackknife resampling method, after 
splitting the UKIDSS DXS area into 48 sub-fields. 

Although our data cover a wide area, it is hard 
to avoid the effect caused by the finite survey 
area, which is referred to as the integral constraint 
(IC, Groth & Peebles 1977) which is additive, i.e., 
w(0) = wiobs (9) + IC. This effect 
rected by the empirical method in 
(Il999h with an assumption of the actual shape 
of two-point correlation function. We apply two 


can be _cor- 
Roche et al. 


different approaches for stellar mass limited sam¬ 
ples in § H] and sSFR binned samples in § [5j For 
stellar mass limited samples, we use each two- 
point correlation function obtained from the halo 
model with observed parameters and then cal¬ 
culate the i ntegral constra i nt wi th the empirical 
equation in Roche et al. I ( 1999h . This process 
is included in the model fitting procedure to ob¬ 
served correlation functions, and taken into ac¬ 
count to find the best fit HOD parameters (see 
Wake et al. 2011). However, the integral con¬ 
straint must be measured differently from stellar 
mass limited samples for sSFR binned samples, 
since we do not perform the halo model i ng for this 
selection. As introduced in Kim et aJ~l ( 2011 1. we 
assume the empirical functional form as the true 
correlation function, and then use the iterative 
method in Roche et al. (199^). In this work, we 


assume the fu nctional form of w (0) = a\9 & 1 + 
o /-2 exp(—/^ 2 ^) ( Kim et al. 20141) . The combina¬ 


tion of a power-law and an exponential compo¬ 
nent is necessary to describe small scale clustering 
of galaxies (power-law component) and large scale 
halo-to-halo clustering (exponential component), 
simultaneously. We note that the integral con- 



























































straints of sub-samples range from 0.001 to 0.013 
and that a more clustered sample tends to have 
a larger integral constraint value. Furthermore, 
for the stellar mass limited samples, the IC values 
measured by using the halo model and the func¬ 
tional form show the differences of <0.0017. 


3.3. HOD Fit 


Since galaxies reside in dark matter haloes, 
the distribution of dark matter haloes or density 
peaks are imprinted on the clustering of galax¬ 
ies. In this context, we are able to link galaxies 
with dark matter haloes through the halo model 
(see Cooray & Sheth 2002 for a review). Here, 
we apply a halo model using the Halo Occupation 
Distribution (HOD) to study the relation between 
galaxies and their host dark matter haloes. This 
model has be en widely app li ed to v arious galaxy 


populations (Zheng et al. 

2007; 

Blake et al. 

2008: Wake et al. 2008| Ross & Brunner 2009a: 

Sawangwit et al. 201 

i: Wake et al. 

2011; Zehavi et al. 

20111 Krause et al. 

20li 

Nikoloudakis et al. 

2013c Kim et al. 20( 

1). 


In order to model the best fit correlation func¬ 
tion, we have to parameterize the basic compo¬ 
nents of the halo model such as the number den¬ 
sity of haloes, the satellite distribution in haloes, 
the halo bias and the mean number of galaxies 
at a given halo mass. First, we adopt the halo 
mass function (n(Mh a in)) and the halo bi as func¬ 
tion (H(Afhaio)) from [Tinker et al.~l ( 201Cll ) for the 
number density and bias of haloes, respectively. 
Second, the distribution of satellites is assumed 
to follow the NFW profile ( Navarro et al. Ifl997t) 
with the concentration parameter depending on 
redshift dBul lock et al. 2001^_ Zehavi_etj3 i L 2004 : 


iBlake et, al. ll2008URoss fc Brunner 112009a ). Since 

ro(Mhaio), B(Mh a io) and the NFW profile depend 
on redshift, we determine these at the mean red- 
shift of our sample, i.e., z = l. Finally, we spec¬ 
ify the mean number of galaxies at a given halo 
mass (N(M\ la \ 0 )), which in turn is parameterized 
for central galaxies ( N c ) and sat elli tes (N,) se pa- 


rately and was introduced bv lzheng et al. (I2005f i , 


as 


N(M halo ) = N c (M ha io) + N s (M hal 0 ), (2) 


with 


N c (M halo ) = 0.5 


1 + erf 


f log 10 (A/halo/Mcut) 

\ & cut 


and 

Ns(Mhalo) — 0*5 


(3) 


1 + erf 


( loglO (Afhalo/A/cut: ) 


( Mialo — 

J V Mo 

(4) 

where Mhaio is a dark matter halo mass. M cut and 
a cu t define the transition halo mass and shape of 
HODs for central galaxies. Mi is the truncation 
mass for satellites, and Mo and a are the thresh¬ 
old halo mass and the slop e for HODs o f satel ¬ 
lites , respectively (see also IZheng et al. (j2005lf 
and Wake et al. ( 201l|P . 

Our survey area of ~140/i _1 Mpc at z = 1 on 
one side is not wide enough to fully constrain all 
five free parameters ( M cu t , a C ut, M p, Mi an d a 


simu l taneously. In previous work (Zheng et al. 


2007 : Brown et al. 1 120081 : IZehavi et al. 2011 ). it 
was reported that Mi is poorly constrained but 
similar to M cut based on the SDSS data. There¬ 
fore we follow the relation, (Mi = M cut ), in 
this work. In addition, M cut can be determined 
by matching the number densit y of galaxies wi th 


given parameters as applied in Ross & Brunner 
(I2009ah . Consequently, we have just three free pa¬ 


rameters {a cu t , Mo and a) to model the real-space 
correlation function. 

We assume that the mean number of cen¬ 
tral galaxies is unity beyond a specific halo 
mass. Recent work based on cosmological sim¬ 
ulations has argued that this can be below unity 
even in very massive haloes due to active g alac¬ 
tic nucleus (AGN) feedback (Kim et al. _ 2009; 


Gonzalez-Perez et al. 12011 : Contreras et al. 12013 1 
if the stellar mass threshold is chosen to be very 
large. However, the effect of AGN feedback is still 
controversial and there is a debate as to whether 
this is positive or negative on short timescales 
(< 100 Myr). Therefore we do n ot consider this 


effect in this work. Furthermore, Zentner et al. 


( 2013h pointed out that any halo model that does 
not take into account the assembly bias leads to 
a systematic error on the fitted result. Although 
this effect is significant for red populations, it is 
much reduced when considering all galaxies. 

Using the basic com ponents ment ione d above, 
we follow the scheme of [Ross fc Brunner" (2009a). 
Briefly, we model power spectra contributed by 
galaxies in the same halo (1-halo term) and in 


Mi 
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different haloes (2-halo term). The 1-halo term 
is distinguished into central-satellite and satellite- 
satellite pairs again. In order to consider the 
underlying dark matter, the matter power spec¬ 
trum is generated by th e ‘CAMB’ software pack¬ 
age dLewis et al. 200(J) including the formulae of 


Smith et ah 02003 ) to model nonlinear growth. 


Then the modeled power spectrum is transformed 
to the real-space corr elation function using the 


Limber equation (Limber 19541 ) to project the 


modeled real-space correlation function to angu¬ 
lar space. The redshift distribution of each sub¬ 
sample is generated not by the best fit photomet¬ 
ric redshift, but by the possible redshift distribu¬ 
tion of each galaxy. We adopt 90% redshift ranges 
from the Le Phare software for each galaxy, then 
assume the Gaussian distribution above and below 
the redshift separately. This distribution of each 
galaxy is stacked to produce the redshift distri¬ 
bution. Finally, the modeled angular correlation 
function is fitted to the observed correlation func¬ 
tion with the covariance matrix derived in § 13.21 to 
find the best fit parameters. We fit the modeled 
correlation function to the observed one over the 
angular range 3"2 < 9 < 0°.7, where the influence 
of the integral constraint is minimized. 

Additionally, the effective halo mass ( M e ff ), 
the effective bias ( b g ) and the fraction of central 
galaxies (f cen ) are derived from the best fit pa¬ 
rameters with 

Meff = J dMhaloMhalon(Mhalo)fV(Mhalo)/flg, (5) 

b s = j dMhaioB(Mhaio)n(Mhaio)A r (Mhaio)/n 9 (6) 

and 


/cen — / dM h alon(A/halo)A r c(A4alo)/n g 


(7) 


where n g is the number density of galaxies, which 
is fixed to the observed value in this work. 

We perform the halo modeling for only galax¬ 
ies selected above a stellar mass threshold, since 
this model is appropriate for mass or luminosity 
limited samples. For the galaxies in stellar mass 
bins, the difference of HODs between two mass 
thresholds is calculated (see § 14.21) . 


4. CLUSTERING WITH STELLAR MASS 

In this section, we describe the results from 
fitting the halo model to galaxies split by dif¬ 


ferent stellar mass thresholds of M* > 10 lo °, 
10 10 ' 5 and 1O 11O M 0 IS 14.11) . Then the HODs for 
galaxies in different stellar mass bins (lO loo M 0 < 
M* < 10 lo - 5 A7 Q , 1O 1O5 M 0 < M* < 10 110 A7q and 
M* > 10 ll o A7g) are compared in § 14.21 Finally 
based on these HODs, the ratio between stellar 
mass and halo mass is discussed in § 14.31 


4.1. Mass Threshold 


In order to investigate the dependence of the 
HOD parameters on stellar mass, we split the se¬ 
lected galaxies at stellar mass thresholds of M* > 
10 10 , 10 10 ' 5 and lO n M 0 . In total, 66,864, 29,250 
and 4,564 galaxies are selected for each stellar 
mass threshold, respectively. The upper panels of 
Figure[G]show the measured angular two-point cor¬ 
relation function for each sub-sample (points with 
error bars). The relatively large error for M > 
1O 11O M 0 comes from the significantly smaller 
number of galaxies. The solid line is the best fit 
halo model. As seen in Figure [Gl the halo model 
reproduces the observed angular correlation func¬ 
tion well. It also shows a relatively good fit even 
beyond the fitted range (0 > 0°.7) and a more pro¬ 
nounced break between the 1- and 2-halo compo¬ 
nents as the balance of central to satellite galaxies 
changes with stellar mass. The best fit parameters 
and derived quantities are listed in Table CD 

The lower panels of Figure [5] display the best fit 
HODs for each sub-sample (solid line). The dotted 
and the dashed lines are HODs for central galaxies 
and satellites, respectively. It is clear that HODs 
for lower stellar mass galaxies extend down to a 
lower halo mass regime. This is well described by 
the best fit parameters of two halo mass thresh¬ 
olds ( M cut and Mo), which decrease for low mass 
galaxies. In addition, the HOD of more massive 
central galaxies shows a gentler transition shape 
compared to less massive centr a l gala xies as al¬ 
ready reported bv IZehavi et al. I ( 2011 ). The best 
fit parameter <r cut describing the transition shape 
is 0.5 for M, > 1O 1O M 0 and 0.6 for M* > 1O U M 0 . 
In order to check the effect of the galaxy number 
density to the best fit parameters, we perform the 
fit with four free parameters ( a cu t , Mo, a and Mi). 
In this case, n g does not directly constrain Mi, but 
is calculated by the best fit parameters for match¬ 
ing the clustering only. Although, the parameters 
show different values from 3 parameter fits, the 
trend is exactly same. These results are presented 
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Fig. 6.- The upper panels show angular two-point correlation functions of galaxies selected by stellar mass 
thresholds indicated in each panel (points with error bars). The solid line is the best fit halo model. The 
lower panels present the best fit HODs (solid line). Dotted and dashed lines indicate HODs for central 
galaxies and satellites, respectively. It is clear that more massive galaxies are in more massive haloes on 
average. 


in Appendix [A] 

Figure [7] shows derived quantities as a function 
of stellar mass thresholds with filled circles. The 
top and middle panels show the effective halo mass 
and the bias for each sub-sample. Clearly massive 
galaxies tend to reside in massive haloes with high 
bias values, consistent with previous work. From 
the bottom panel of Figure [3 we can see that mas¬ 
sive galaxies also tend to be the central galaxy in 
a massive halo. 


Recently, many researchers have measured the 
clusteri ng of galaxies sp l it by stellar mass. For in¬ 
stance, (hbucaudetj^ ( 20101 1 measured the clus¬ 
tering of galaxies from the Palomar Observatory 
Wide-held Infrared Survey (Conselice et al. 2007). 
Their results show a halo mass of ~ 10 12 S h~ 1 M <3 
and a bias of 1.9 for lO n M 0 < M* < 1O 12 M 0 
galaxies at 0.8 < 2 < 1.2, which is lower than 
ours (cross in Figure 0. However, their analy¬ 
sis was done over a 1.5 deg 2 area, which is >5 
times smaller than ours. In addition, they mea¬ 
sured halo masses based on the “one galaxy per 
halo” assumption with a correction for the halo 


occupation. 

As mentioned above, HOD analysis may be af¬ 
fected by cosmic variance if the survey area is 
too small and the quality of the photometric red- 
shifts used is poor. With these issues in mind we 
quote two more results, one from a much wider 
area and another on e based on the spectro scopic 
information. First, ICounon et al. ( 2012 1 per¬ 
formed the halo modeling to reproduce the clus¬ 
tering of galaxies categorized by luminosity and 
type from the full CFHTLS-Wide survey area 
(~ 133 deg 2 ). Of our sub-samples, galaxies with 
M* > 10 1o ' 5 Mq and 1O 11O M 0 have similar num¬ 
ber densities to their sub-samples containing all 
galaxies at M g — 51ogh < —20.8 and < —21.8 
at 0.8 < 2 < 1.0. Their estimates were M e ff of 
10 12 - 95 ft- 1 M© and 10 13 ' 25 bias of 1.74 and 
2.25, and the satellite fraction (f sa t) of 0.13 and 
0.06 for M g — 51ogh < —20.8 and < —21.8, respec¬ 
tively. Overall, these results are in good agreement 
with our estimates, although the halo mass for the 
brightest luminosity bin is about 0.15 dex higher 
than ours (open squares in Figure 0. Second, 
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Fig. 7. Derived effective halo masses (top), 
biases (middle) and fractions of central galaxies 
(bottom) from the halo modeling. More massive 
galaxies reside in more massive haloes and tend 
to be cent ral galaxies. As r eferen c es, previous re¬ 


sults f rom iFoucaud et al. (20101. Coupon et al. 
( 20121) and Mostek et al. ( 2013i ) are also dis¬ 


played. For the display purpose, literature values 
are slightly shifted on the M* axis. 


Mustek et ah ( 20131) used galaxies with spectro¬ 


scopic inform ation from the DEEP 2 Galaxy Red- 
shift Survey ( Newman et al. EqU to measure a 
projected correlation function of galaxies at 0.74 < 
z < 1.4 at different stellar masses, SFRs and sS- 
FRs. They estimated bias and mean halo mass as 
(1.77, 1O 12 ' 95 /i _ 1 M 0 ) and (1.87, 10 13 - O5 ft- 1 M 0 ) 
for all galaxies at 0.74 < z < 1.05 with M* > 
10 10 ' 5 and 1O 1O 8 M 0 , respectively (open triangles 
in Figure [7]). Although the stellar mass thresholds 
are slightly different from our samples, all esti¬ 
mates are in good agreement with our values. 

It is interesting to note that our estimates 


are consistent with CFHTLS (wide-area survey) 
and DEEP2 (a wide spectroscopic survey), while 
small area su rveys give somewhat different results. 
Kim et al. (12011 ) demonstrated the importance 
of surveyed area for the reliable measurement of 
galaxy clustering on large scales in overcoming 
cosmic variance studying Extremely Red Objects. 
The scatter of clustering on large scales is ~30% 
and ~20% in areas of 0.26 deg 2 and 0.6 deg 2 , re¬ 
spectively. In addition, the halo model with HODs 
may be a more reliable scheme to measure the 
masses of host haloes for all, central and satellite 
galaxies, since this takes into account a more real¬ 
istic contribution of central and satellite galaxies. 
Based on HODs discussed in this section, we can 
extend our analysis for galaxies split into stellar 
mass bins, which is an aspect that has not been 
fully explored in the literature. 


4.2. Mass Bin 

If we are able to find the HODs of galaxies in 
various stellar mass bins, we can understand the 
relation between galaxies and their host haloes. 
However it is difficult to model the clustering of 
galaxies in stellar mass bins, since a different HOD 
shape should be assumed that is not well under¬ 
stood yet. In this section, we use the HODs dis¬ 
cussed in § 14.11 to obtain the HODs of galaxies in 
different stellar mass bins. 

Zehavi et al. I ( 201ll) took the difference of 


HODs of luminosity threshold samples, in or¬ 
der to obtain the HODs of luminosity binned 
samples. We adopt this method to obtain the 
HODs of galaxies in different stellar mass bins. 
For instance, we obtain a HOD of galaxies at 
1O 1O '°M 0 < M* < 1O 1O ' 5 M 0 by subtracting the 
HOD of M* > 1O 1O ' 5 M 0 galaxies from the HOD 
of M* > 1O 1O M 0 galaxies. The bottom panels of 
Figure [H] show the HODs for 1O 1O °M 0 < M* < 
10 1O - 5 M 0 (left) and 1O 1O - 5 M 0 < M* < 1O 11O M 0 
(middle) galaxies, respectively, estimated in this 
way. For comparison, we also plot the HOD for 
M* > 10 11 -°M 0 galaxies in the right panel. The 
upper panels display the observed angular cor¬ 
relation functions (points) for each stellar mass 
binned sample. The solid lines in the top-left 
and top-middle panels are not the best fit re¬ 
sults, but the modeled clustering based on the 
stellar mass binned HODs. The top-right panel is 
identical to that in Figure [6] The modeled cor- 
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Fig. 8.- The similar plot with Figure El for stellar mass binned galaxies. All symbols are the same to 
Figure [6] However, the solid lines in the upper panel are reproduced angular correlation functions from the 
calculated HODs based on those for mass threshold galaxies (see text for more details). The most massive 
bin is identical to the right panel in Figure El More massive central galaxies are more widely distributed in 
halo masses than less massive ones. 


relation functions show a good agreement with 
the observed ones, except at a very small scales 
(6 ~ 0°.001) where the modeled correlation func¬ 
tions are slightly underestimated in comparison to 
the observed ones. 


Returning to the HODs, we note that the HOD 
for central galaxies is similar to a log-normal dis¬ 
tribution, which indicates that there is a crude cor¬ 
relation between the host hal o mass and the stella r 
mass of the central galaxies ( Moster et al. 20ldl . 


The width of the distribution reflects the amount 


of the scatter in this relation (jZheng et, al. 1120051) . 


and we find that the central HOD for more mas¬ 
sive galaxies covers a wider halo mass range at a 
given value of N c (Mh a io ) than that for less mas¬ 
sive ones with the peak position shifted toward 
the high halo mass regime. As shown in the previ¬ 
ous section, more massive galaxies show a gentler 
central HOD shape than less massive ones, which 
means a larger scatter between the stellar mass of 
central galaxies and the host halo mass. Therefore 
the broader shape of the HOD of massive galax¬ 
ies reflects this effect. This can be caused by the 


stellar mass growth being stopped by any quench¬ 
ing mechanism, whilst_the halo es keep grow¬ 


ing. For instance, iGonzalez-Perez et al. ( 201li ) 


pointed out that the inclusion of AGN feedback 
leads to a change in slope and a larger scatter 
of the relation between AT-band luminosities of 
central galaxies and their host halo masses, since 
the feedback presents gas cooling. In addition, 


Contreras ct al. , ( 20131) compared the HODs pre¬ 


dicted by different semi-analytic models and show 
that the central HODs of more massive galaxies 
are more affected by this feedback in all simula¬ 
tions. The idea is that above a certain luminosity 
or mass, galaxies do not grow to have a very large 
stellar mass due to the negative feedback by AGN, 
resulting in cases where very massive haloes pos¬ 
sess central galaxies with reduced (but still mas¬ 
sive) stellar masses. Consequently, the HOD of 
massive galaxies becomes extended toward a large 
halo mass. 


The upper panel of Figure El shows the effec¬ 
tive mass of dark matter haloes hosting galaxies in 
each stellar mass bin calculated from the HODs. 
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Fig. 9.- The upper panel shows the derived 
halo masses for central (filled circles) and satel¬ 
lite galaxies (open circles) in different stellar mass 
bins. The bottom panel represents the fraction of 
central galaxies. For the highest stellar mass bin, 
we use the M* > lO n M 0 galaxies. The horizontal 
bars indicate the size of stellar mass bins. 


For the most massive bin, we plot the stellar mass 
threshold sample of M* > 10 11 Mq. The mass of 
dark matter haloes hosting central galaxies (filled 
circles) increases as the stellar mass increases with 
close to a linear slope. That for satellites (open 
circles) also shows a similar trend, but not as dra¬ 
matic as it is for central galaxies. The reason for 
this difference in satellite galaxies is that a mas¬ 
sive halo contains a large number of satellites in 
addition to a massive central galaxy. The bottom 
panel of Figure [9] displays the fraction of central 
galaxies in each stellar mass bin. This indicates 
that the massive galaxies are more likely to be 
central galaxies as the halo mass increases. 

4.3. Stellar Mass to Halo Mass Ratio 

Since galaxies evolve in dark matter haloes, 
the properties of galaxies depend on their host 
dark matter haloes. In this context, the rela¬ 
tion between stellar mass (M*) and halo mass 


(T/haio) is a good testbed to constrain the evo¬ 
lution of galaxies. In addition, the ratio be¬ 
tween the stellar mass of central galaxies and 
the halo mass is sensitive to the conversion effi¬ 
ciency from baryons to stellar m ass in the cen¬ 
tral galaxy ( Zheng et al. I 120071) . In previous 
work, the M*-M ha i Q (or M*/M halo -M halo ) rela¬ 
tion was derived using several different methods 
such as HOD, Conditional Luminosity Function 
and Sub-halo Abundance Matching rtZheng et al. 
2007 : Behroozi et al. 2010t Moster et al. 2010t 


Wang fc .Ting ll201flT Wake et al. 1 201 lh . However. 

one of the advantages of the HOD framework is 
that the clustering and the number density of 
galaxies are fitted simultaneously without any as¬ 
sumption about the scatter between M* and Mh a i 0 
being a fixed value. Here, we derive the relation 
based on the best fit HODs directly, unlike previ¬ 
ous HOD work, fitting a functional form to best 
fit HOD parameters. 


We use the HODs for stellar mass binned sam¬ 
ples discussed in the previous section to calculate 
the stellar mass to halo mass ratio (M*/Mh a i 0 ). 
Since the HOD is the mean number of galaxies at 
a given halo mass, the summation of HODs mul¬ 
tiplied by mean stellar masses for a given stellar 
mass bin represents the stellar mass at a given halo 
mass. Therefore, the stellar mass at a given halo 
mass can be calculated by 


M*(M halo ) = ^ JVi(M halo )(M.)i ( 8 ) 


where iVj(Mh a i 0 ) is the central or satellite HOD for 
the ith stellar mass binned sample and (Af*)j is the 
mean stellar mass of galaxies in the ith stellar mass 
bin. In this work, there are three stellar mass bins, 
each corresponding to 1O 1O M 0 < M* < 1O 1O ' 5 M 0 , 
lO ia5 M 0 < M* < 1O 11O M 0 and M* > 10 n M e . 
Using this stellar mass, we calculate the stellar 
mass to halo mass ratio. 


The upper panel in Figure [TO] shows the stel¬ 
lar mass to halo mass ratio for the central galaxy. 
The red solid line is the result based on equa¬ 
tion (8). The peak of the ratio is located at 
~ 10 1216 /i _1 M 0 with a ratio of 0.024. Below and 
above this halo mass, the ratio drops rapidly. This 
means that the conversion from baryons to stellar 
mass in the central galaxy is the most efficient 
in ~ lO 12 h -1 M 0 dark matter haloes as traced at 
z ~ 1. For the comparison, we also disp la y the 


ratio at 2 = 1 (royal blue line) from Moster et ah 
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( 2010h . We have to note that the ratio may be 
underestimated for low mass haloes in this work, 
since the stellar mass limit is M* = 10 1o Mq, and 
less massive galaxies are missed. Therefore the ac¬ 
tual slope of the ratio in the low halo mass regime 
may be flatter than our measurement. 



Fig. 10. The stellar mass to halo mass ratio 
based on the best fit HODs. The upper panel 
is for central galaxies only at z ~ 1. The re¬ 
sult of this work is represented by a thick red 
line. For the compariso n, we also plot previ¬ 
ous results at z ~ 1 from Behroozi et al. I ( 201dh 
(green clotted line), Foucaud et ah fcoidi (black 


points), iMoster etak 
line) and]Wake et ah 


1 20101) (blue short dashed 
(l201ll) (cyan points). The 
ratio shows a peak at logMhaio ~ 10 12 h _1 M Q , 
which is consistent with previous results. The 
orange line is the ratio calculated by the empir¬ 
ical function in lWake et al. ( 2011 ) for our result. 
The lower panel displays the ratio for central (red) 
and satellite (blue) galaxies. The black line is the 
combination of both central and satellite galax¬ 
ies. The orange and cyan lines are results with 
varying HOD parameters randomly within the 
uncertainty range for central and satellite galax- 
ies, respective l y. So lid and dashed lines are from 
Moster et al. (2010) for 2 =0 and 2.5 respectively. 
Satellites mainly contribute to the total stellar 
mass in group or cluster scale haloes. 


Here, we al so compare our resu lt with other pre¬ 
vious works. I Wake et al. ( 2011 ) also performed 
the halo modeling for stellar mass limited sam- 


(Brammer et al. 

2009| van Dokkum ct al. 2009. 

2010t Whitaker et al. 201ll). After fitting the ern- 

pirical function in Zehavi et al. ( 

201l|) to the best 


fit HOD parameters, they obtained the ratio at 
z = 1.1 (cyan points in Figure H0l) . Their result 
shows a peak at a higher halo mass, and a flatter 
(steeper) shape in a high (low) halo mass regime 
than ours, although the relation is poorly con¬ 
strained at high and low halo masses. If we also fit 
the same function to our best fit parameters, this 
also shows a steeper slope in the low halo mass 
regime, but the discrepancy is not so significant 
in the high mass regime (orange line). On the one 
hand, this means that the functional form may un¬ 
derestimate the ratio for less massive haloes, since 
this does not fully take into account HODs for stel¬ 
lar mass binned samples. On the other hand, for 
the high halo mass regime, the flatter trend may 
be caused by t he differ ent parameter set for HODs 


as I Wake et al. ( 2011 ) fixed the transition shape 


of central HODs with a cu t = 0.15. In fact, the re¬ 
lation between derived halo mass parameters for 
central gala xies (M r/n t i n equa tion (3)) and stellar 


masses in lWake et al. ( 2011 ) is flatter than ours. 


If we make the relation for our fitted parameters 
flatter arbitrarily, the ratio becomes much fiat- 


ter in the high halo mass regime. iFoucaud et al 


( 2010h also estimated the ratio as 0.032±0.022 and 
0.026T0.017 for lO n M 0 < M* < lO n ' 5 M 0 and 
1O 115 M 0 < M* < 1O 12O M 0 at 0.8 < z < 1.2, 
respectively (black points). Although these are in 
agreement within the uncertainty range, their halo 
masses are smaller compared to ours as mentioned 
above. Therefore it is possible that their points 
move toward a lower ratio and a higher halo mass 
regime, and become consistent with ours. 

Additionally, iBehroozi et al. ( 2010l) obtained 
the same quantity using the abundance match¬ 
ing technique. We plot their result for the cen¬ 
tral galaxy including full uncertainties at z = 1 in 
the upper panel of Figure [TO] (green points and 
line). Our result agrees with theirs within the 
uncertainty. However, our result may be slightly 
steeper i n both the high a n d low halo mass regimes 


than the iBehroozi et ah 


similar discrepancy from IMoster et al. 


(120101) resu lt. We find a 
(120101) as 


well (royal blue). As mentioned above, the ratio 
in the low halo mass regime is underestimated due 
to the stellar mass limit in this work. However, 
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in the high halo mass regime, that is not the case 
and ot her factors may a c count for th is discrepancy. 
First, iBehroozi et al. ( 2010h and Moster et, al 


( 20101 ) included the scatter in the M*-Mi ia io rela¬ 
tion with a fixed value. Whereas, our HOD takes 
into account such a sc atter as shown in Figure [ 6 ] 
Behroozi et al. ( 2010 ) pointed out that a larger 


scatter in the relation makes the relation steeper. 
Second, the stellar mass uncertainty on individual 
galaxies can also affect the result, since the num¬ 
ber of low mass galaxies having overestimated stel¬ 
lar masses is larger than that of high mass galaxies 
with underestimated stellar masses. If this scat- 
ter is included, the rel ation also becomes steeper 
( Behroozi et al. 201dh . 

In addition to the central galaxy, satellites also 
contribu te to the total bary ons i n dark matter 
haloes . Behroozi et al. ( 2010h and lFoucaud et, al~ 


( 2010 h pointed out that satellites account for the 
majority of the total stellar mass in more massive 
dark matter haloes. As we discussed in the previ¬ 
ous section, satellite galaxies with M* > 1 O 1 O M 0 
tend to live in massive haloes. The bottom panel 
in Figure [10] also shows this trend. The red line 
is the ratio for the central galaxy only, which is 
the same as in the upper panel. The blue line 
is for satellites calculated by equation ( 8 ) with 
the HODs in Figured] However, in this case, the 
value represents the ratio of the total stellar mass 
in satellites to the total mass of the dark matter 
halo. Orange and cyan lines are the results for 
central and satellite galaxies, respectively, which 
are derived from 100 different HOD sets with vary¬ 
ing parameters and a ) with values var¬ 

ied randomly within their measured uncertainty 
ranges. In contrast to the central galaxy, the 
ratio for satellites increases with increasing halo 
mass. In addition, the M*/Mh a io ratio for satel¬ 
lites is comparable with the value for the cen¬ 
tral galaxy in group environments (10 13 /i^ 1 Mq < 
Mhaio < 10 14 /i _1 Mq) and dominant in cluster en¬ 
vir onments (M^aio > 10 1A h ~ 1 M c ^ ). As discussed 
Behroozi et al. I ( 2010h and Foucaud et al. 


(1201 oh . the accretion of dark matter toward a 
massive halo leads to the rapid growth of the halo 
mass. At the same time, infalling galaxies become 
satellites and contribute to the total stellar mass. 
However, the growth of central galaxies is not so 
efficient in this regime due to the deep potential 
well preventing an efficient merger or cold gas ac¬ 


cretion onto the central galaxy. We also note that 
our data do not include faint, possibly numerous 
satellite galaxies. Thus, the M*/Mh a i 0 ratio for 
satellites may be underestimated and can only be 
addressed with deeper data. 


From the M*/Mh a io ratio for the central galaxy 
and satellites, we find that the star formation ef¬ 
ficiency for the central galaxy is the most efficient 
in ~ 10 12 1i~ 1 Mq dark matter haloes at z = 1 
with the M*/Mhai 0 ratio peaking at about 0 . 02 . 
Moreover, we show that satellites are the domi¬ 
nant contributor to the M*/Mh a io of high mass 
haloes (Mh a i 0 > lO 13 - 5 h _ 1 M 0 ). When we con¬ 
sider both central and satellite galaxy samples 
with M* > 10 1o Mq, we find the star formation 
efficiency in terms of total halo mass and stars in 
all galaxies is l%-2% consistently. However, the 
cut-off at low halo mass is created by our selec¬ 
tion of massive galaxies, so the mass of stars is 
not fully accounted for in the lowest mass haloes. 
We also confirm the evolution of this relation, 
comparing to t he res ult at different redshift from 
Moster et al. ( 201 0) (the bottom panel in Fig¬ 


ure [lO]). 


5. CLUSTERING WITH STAR-FORMING 
ACTIVITY 

The clustering of galaxies has been shown to 
depend on star formation rate and color. In this 
section, we split galaxies at 0.8 < 2 < 1.2 based 
on their star forming activity to investigate how 
their clustering properties depend on star forma¬ 
tion. Additionally, since sSFR is roughly related 
to the star formation efficiency in a galaxy, we use 
various sSFR criteria to define each sub-sample. 

5.1. Passive vs. Star Forming 

We start with the comparison of bias values 
for different populations of galaxies. For this 
analysis, galaxies are split into passive and star¬ 
forming. In order to define the sSFR cut for pas¬ 
sive galaxies, we check galaxies detected by Spitzer 
from the IRSA catalog. We cross-match Spitzer 
sources to our galaxies with M* > 10 lo M Q and 
0.8 < z < 1.2, and find 1871 IRAC and 213 
24/um sources. Of 24/zm sources, 15% and 8 % 
satisfy logsSFR/yr _1 < —10 and —10.5, respec- 


3 http: //irsa. ipac.caltech.edu/applications/Gator/ 
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tively. This means that the sSFR cut may not 
be a clean way to classify pure passive galaxies, 
but we can isolate them efficiently. Thus, al¬ 
though there are two galaxy loci in Figure [5] sep¬ 
arated by logsSFR/yr _1 ~ —10, we classify pas¬ 
sive galaxies with logsSFR/yr _1 < —10.5 in order 
to minimize the contamination by dusty galax¬ 
ies. Then, star-forming galaxies are defined by 
logsSFR/yr” 1 > —10. we then apply stellar mass 
cuts with M, > 10 10 0 , 10 10 ' 5 and 1O 11O M 0 for 
each population. However, since the number of 
star-forming galaxies in the highest mass bin is 
too small to measure the clustering, we do not in¬ 
clude this sub-sample in our analysis. 
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Fig. 11. Biases for passive (red) and star¬ 
forming (blue) galaxies in logM*/M 0 > 10.0, 10.5 
and 11.0 (filled circle s). Open triangles are from 
[Hartley et al. I ( 2013ll with 10.0 < logM*/M 0 < 
10.5 and logM»/M 0 > 10.5. Open squar es are 
for red galaxies in ICounon et al. I ( 2012h with 
M g — 51ogh < —20.8 and —21.8. Also the open 
circle shows the result_of LRGs at z ~ 1 
Nikoloudakis et al~l ( 20131) . 


z ~ i m 
For the display pur¬ 
pose, open symbols are slightly shifted along the 
M* axis. 


For each selected sub-sample, the angular two- 
point correlation function is measured. Since the 
intrinsic HOD shape for each sample is not well 
understood, we estimate the bias by fitting the an¬ 
gular correlation function of the underlying dark 
matter instead of finding the best fit HODs. First, 


we obtain the real space correlation function of 
dark matter from the linear matter power spec¬ 
trum at z = 1. Then the real space correla¬ 
tion function is transformed into the angular cor¬ 
relation function with the observed redshift dis¬ 
tribution, after multiplying it by the bias value. 
Through fitting the transformed angular correla¬ 
tion function to the observed one, we find the best 
fit bias values for each sub-sample. The fitting is 
performed to the angular range between 0.02° and 
0.7°, where the influence by the integral constraint 
is minimized and 2-halo component dominates. 


Filled circles in Figure Qj] show the estimated 
bias for passive (red) and star-forming (blue) 
galaxies. It is clear that passive galaxies are 
more clustered than star-forming ones, which 
is consistent with pre vi ous r esult s at the si mi¬ 
lar redshift ( Coil^aL _ 200t|_ McCrackenet_al 


2008: 1 Williams et al. 120091 : lHartlev et al. Il2010l 
2013 : lBielbv et al 1120141) . Additionally, high mass 


galaxies show a stronger clustering strength than 
low mass ones independently of population. On 
the other hand, some previous results pointed 
out that passive or red galaxies show similar (or 
higher) clustering strengths with decreasing stellar 
masses (DEEP2 in Coil et al. 2008 and UKIDSS 
UDS in Williams et al. 2009 and Hartley et 
al. 2013). Open triangles in Figure |TT] are for 
galaxies with 1O 1O M 0 < M* < 10 1 05 Ai^ an d 
M* > lO lo ' 5 Af 0 at z ~ 1.06 from Hartley et al 


( 20131) based on the UKIDSS UDS data. Pas¬ 
sive galaxies show a similar clustering strength in 
both bins, which implies that there is no stellar 
mass dependence. This was interpreted as the 
contribution of low mass satellite galaxies in mas¬ 
sive haloes, leading to similar clustering strengths 
independently of the stellar mass. 


Our results may appear to be different from pre¬ 
vious literature results. However, we also display 
bias values for red galaxies with M g — 51ogh < 
—20.8 and —21.8 (re d open square) a t 0.8 < z < 
1.0 from CFHTLS in I Coupon et al.l ( 20121) . The 
number densities in each magnitude bin are similar 
to those for our passive galaxies with M* > 10 10 5 
and 10 11 Mq, respectively. Although their mea¬ 
surements are lower than ours, brighter red galax¬ 
ies are more clustered than fainter ones, which 
is the same trend to our passive galaxies. Pos¬ 
sibly lower values are caused by different selection 
criteria, since absolute magnitude limited samples 
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Fig. 12.— Angular correlation functions (points) for sSFR binned samples. The top row is for passive 
galaxies (logsSFR/yr -1 < —10.5) with different stellar mass thresholds. The middle and bottom rows are 
for star-forming galaxies in various sSFR bins as labeled. The solid line is the best fit result, and the dotted 
line is for M* > 1O 11 M 0 passive galaxies. Comparing solid and dotted lines, the amplitude difference is the 
largest in the bottom-left panel, which indicates the lowest bias. 


can include low stellar mass galaxies, but exclude 
some high mass ones. In addition, they noted 
that red galaxies at z > 0.8 were contaminated by 
bl ue galaxies. The similar trend was also found 


McCracken et al. . ( 2008i ). which showed early- 


type galaxies with brighter absolute magnitude 
are more strongly clustered than fainter early-type 
galaxies at 0.7 < 2 < 1.1, based on the CFHTLS 
data. Finally, we also n ote t he bias value from 
LRGs at z ~ 1 from Nikoloudakis et al. ( 2013 1 


(red open circle). The number density of LRGs in 
their paper is similar to that of our high mass pas¬ 
sive galaxies, and the estimated bias is also con¬ 
sistent. 


Since the area coverage of our data is not wide 
enough to fully overcome the effect of cosmic vari¬ 
ance, it may be hard to conclude whether the clus¬ 
tering strength of passive galaxies at z ~1 is corre¬ 
lated with stellar mass in our dataset conclusively. 
However, our result shows the consistent trend 
with the result from the wider survey data, which 
indicates our dataset is not significantly affected 
by cosmic variance. Based on this conclusion, we 
investigate the clustering property of galaxies with 


finer sSFR bins below. 

5.2. sSFR bins 

As shown in Section 13.1.21 the main sequence 
of star-forming galaxies has a slope of —0.33 on 
the sSFR versus stellar mass plane. However, our 
main goal is to investigate the clustering strength 
as a function of star formation activities at 2 ~ 
1. Thus, we simply apply five sSFR bins for 
star-forming galaxies (—8.5 < logsSFR/yr _1 < 
-8.0, -8.85 < logsSFR/yr" 1 < -8.5, -9.25 < 
logsSFR/yr _1 < —8.85, —9.7 < logsSFR/yr -1 < 
-9.25 and -10.1 < logsSFR/yr" 1 < -9.7). The 
width of each bin is determined to include suf¬ 
ficient galaxies for a reasonable clustering mea¬ 
surement. Additionally, passive galaxies are de¬ 
fined as galaxies with logsSFR/yr -1 < —10.5, 
and they are further divided into mass thresholds 
of M* > 10 10 , 10 10 ' 5 and 10 n M Q . Column 3 in 
Table m lists the number of galaxies in each bin. 
We also compare the clustering properties in each 
sSFR bin with narrower stellar mass bins later. 

Figure [12] displays the angular correlation func- 
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tion of each sub-sample (points) and the best fit 
result (solid line). The top three panels are for 
passive galaxies with different stellar mass thresh¬ 
olds. The second and third rows are results for 
star-forming galaxies in various sSFR bins with 
M* > 1O 1O M0. The labels in the second and 
third rows indicate the sSFR range in a loga¬ 
rithmic scale. We also plot the best fit result 
for M* > 10 11 M q passive galaxies with dotted 
line. After comparing the solid and the dotted 
lines for sSFR binned samples, the difference is 
the largest for —8.85 < logsSFR/yr _1 < —8.5 
galaxies, which means it has the lowest bias among 
our sub-samples. On small scales, we are able 
to find excess of clustering amplitudes from all 
sub-samples, although the measurement with rel¬ 
atively large uncertainties is affected by the small 
number of objects. However, the excess for star¬ 
forming galaxies looks less significant than passive 
ones. This may mean that the number of star¬ 
forming satellites in each bin is lower than passive 
satellites, which intrinsically weakens the cluster¬ 
ing strength. However, the actual spatial distribu¬ 
tion of satellites in each sub-sample may also influ¬ 
ence the clustering. For an extreme example, if a 
star-forming sub-sample includes one satellite per 
halo and just a portion of centrals, the small scale 
clustering is more weakened than the large scale. 
Unfortunately, our data is not enough in depth 
and area to demonstrate these effects separately 
(or simultaneously). Larger and deeper datasets 
in the future will allow us to investigate this issue. 
In this work, we will focus on the bias or halo mass 
estimated from the large scale clustering only. 

Figure [13] shows the estimated bias as a func¬ 
tion of sSFR (points). The size of the points for 
passive galaxies represents the stellar mass thresh¬ 
olds which are M* > 10 10 , 10 10 5 and 1O 11 M 0 
from small to large points, respectively. For 
star-forming galaxies, a single mass threshold of 
M* > 10 1 o Mq is applied. For comparison, we also 
plot the result of blue gala xies with Mb < —20.5 
at 0.74 < 2 < 1.05 from IMostek et al. (f2013h 
(open squares). As found by IMostek et al. 

'* 201:1 ). we also confirm that the bias or a clus¬ 
tering strength decreases with increasing sSFRs 
up to logsSFR/yr _1 ~ —8.6. The discrepancy 
of biases may come from different selection cri¬ 
teria, since they have used rest-frame magnitude 
limited samples that may include M < 10 lo Mg 



Fig. 13.- The bias of sSFR binned galaxies. The 
vertical dotted line shows the criterion distinguish¬ 
ing star-forming and passive galaxies. For passive 
galaxies (logsSFR/yr^ 1 < —10.5), a symbol size 
corresponds to stellar mass thresholds, M* > 10 10 , 
lO 10 5 and 10 n Mg from small to la rge p oints. The 
previous results for z ~ 2 sBzKs dLin et al. 2012 1 


and 2 ~ 1 galaxies (jMostek et ah 120131 1 are also 


plotted with open triangles and open squares, re¬ 
spectively. The clustering strengths increase with 
decreasing sSFRs. However, it is also found that 
the highest sSFR galaxies are more strongly clus¬ 
tered than main sequence galaxies. 


galaxies which can dilute their measured clus¬ 
tering strength. Since the sSFR value for the 
main sequence at 2 = 1 is logsSFR/yr _1 ~ —9 
in Figure [5] galaxies just above (or the upper 



sible reversal of the relation in the highest sSFR 
bin (—8.5 < lo gsSF R/yr ^ < — 8.0). This was al¬ 
ready noted bv lLin et al. ( 2012 ) for sBzK galaxies 
at 2 ~ 2, but we find a similar trend at 2 ~ 1 and 
we plot their results with open triangles. Direct 
comparison of the absolute values is difficult due to 
different stellar mass limits. Galaxies in this work 


are more massive than those in [Lin et al. ( 201211 . 


However, both results show the same trend, and 
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Fig. 14.— Angular correlation functions for sSFR binned samples in 10 10 < M*/M 0 < 10 10 ' 5 (blue) and 
lO 10 5 < M*/M 0 < 10 11 (red). Points with error bars are correlation functions measured and solid lines are 
the best fit dark matter clustering on large scales. Black labels in each panel indicate the sSFR range in a 
logarithmic scale. 


we find a similar sSFR for the lowest bias sample. 

Additionally, we note the mass of haloes hosting 
each sub-sample. The halo mass is calculated us- 
the measured bias and the halo mass function 
in lTinker et al. I ( 20ld b Since this procedure does 
not consider the inclusion of satellites, the calcu¬ 
lated halo mass may show a discrepancy from that 
from the full HOD framework in previous sections. 
However, since our main goal here is to examine 
how the bias factor or host halo mass change as a 
function of sSFR, we will defer an HOD analysis to 
a future paper where we will include more galaxies 
over a wider area. Details of the influen c e of s atel- 
lites are well described in Zheng et al. ( 2007h . In 
Appendix [5] we also compare halo masses derived 
by the halo model and the direct fit. Our mea¬ 
sured biases and halo masses are given in Table [21 
We find that the mass of host haloes of galaxies 
with the lowest bias is Mhaio = 10 12 - 684 /i _1 M 0 . 
Also the mean halo masses for passive galaxies 
are M halo = 10 12817 , 10 12 ' 940 and 10 13207 h,- 1 M Q 
from the lowest to the highest stellar mass bins, 
respectively. 

From Figure fldl we show the clustering strength 
as a function of sSFRs. Now we investigate the 


same quantities using galaxies in different stellar 
mass bins, in order to check the influence of stel¬ 
lar mass as the negative main sequence slope on 
the sSFR-mass plane could arise from a selection 
bias toward high sSFR galaxies that are in less 
massive haloes and hence are less clustered. To 
assess the level of this potential bias we further 
split galaxies into different stellar mass bins of 
1O 1O M0 < M* < 10 1o " 5 Mq (low mass, LM) and 
10 1O 5 Mq < M* < lO n M 0 (high mass, HM) with 
the same sSFR bins mentioned above. Unfortu¬ 
nately, it is hard to measure the secure angular 
correlation function of star-forming galaxies with 
M* > 1O 11 M 0 divided into various sSFR bins be¬ 
cause of the small number of objects. Thus, we 
do not consider the highest stellar mass bin here. 
Figure IT4l shows the observed correlation function 
of sSFR binned samples in LM (blue) and HM 
(red). Symbols are the same as Figure [T2l Black 
labels indicate sSFR ranges. As already shown 
in §EU passive galaxies in HM is more clustered 
than those in LM. On the other hand, star-forming 
galaxies with log sSFR < —9 in LM show stronger 
clustering strengths than HM ones, but this is 
opposite at the high sSFR range. 

Figure [15] displays the halo mass as a function 
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of sSFR for the 1O 1O M 0 < M* < 1O 1O 5 M 0 (blue) 
and lO la5 M 0 < M* < lO n M 0 (red) bins, respec¬ 
tively. In this case, we plot the halo mass on the y- 
axis rather than the bias factor. From both stellar 
mass bins, we confirm the anti-correlation in the 
sSFR-Mhaio relation in the low sSFR regime of 
star-forming galaxies and at lower significance the 
reversal of the relation at logsSFR/yr _1 > —8.5 
shown in Figure Q21 This means that the feature 
shown in Figure [13] is mainly caused by an envi¬ 
ronmental effect connecting the halo occupation of 
high sSFR galaxies to the mass of their host halo 
and not due to an effect intrinsic to their stel¬ 
lar mass. One difference between two mass bins 
is the sSFR value showing the lowest halo mass. 
Although the uncertainty is substantial, HM star¬ 
forming galaxies with the lowest halo mass have 
lower sSFR than LM ones. Considering the slope 
for main sequence galaxies in Figure [5] galaxies on 
the upper envelope of the main sequence in each 
stellar mass bin reside in the lowest mass haloes. 
The number of galaxies in each bin and all the 
measured values are listed in Table [3] 


Finally, it is worth comparing our results to 
the known sSFR-lo cal density relation. First, 


Cooper et ah (120081) found the same anti-correlation 


at z ~ 1 over the similar sSFR range from the 
DEEP2 survey. In addition, we can also see a sug¬ 
gestion of a weak reverse relation in the highest 
sSFR bin ( see F igure 9 in their paper). Second, 
Ko et al. ( 2012h also reported the similar trend 


at z = 0.087 based on the AKARI North Eclip¬ 
tic Pole (NEP) Wide survey and A2255 from the 
AKARI CLusters of galaxies EVoLution studies 
(CLEV L; Im et a l. 200 8; Lee et al. 2009a). Al¬ 
though IKq et al~ ( 20121 1 focused on the very low 
redshift, there is no significant difference in the 
bias-sSFR relations at low and high redsliifts as 
we already mentioned above. Moreover, their field 
coverage is wide enough to study the relation from 
low to high local densities. In conclusion, our re¬ 
sults also show good agreement with local density 
studies. 


6. DISCUSSION 

In the previous section, we found a clear de¬ 
pendence of clustering strength on sSFR. Here, we 
further investigate the property of galaxies in the 
highest sSFR bin, which are in denser environ- 



log sSFR (yr _1 ) 


Fig. 15. The similar plot with Figure [13] for 
1O 1O M 0 < M* < 1O 1O ' 5 M 0 (blue open circle) 
and 1O 1O5 M 0 < M* < lO n M 0 (red filled cir¬ 
cle) galaxies. In this plot, we quote calculated 
halo masses instead of bias values. We confirm 
the trend in Figure 1131 from narrower stellar mass 
binned samples. The arrows present possible evo¬ 
lutionary paths. See text for more details and 
discussion. The horizontal dashed lines indicate 
M eff values of haloes for all galaxies in each stel¬ 
lar mass bin. Downward arrows at the bottom are 
central sSFR values of the main sequence for each 
stellar mass bin. 


ments than main sequence star-forming galaxies, 
and we present a possible scenario for the evolu¬ 
tion of star-forming galaxies. 

We begin with the comparison of the internal 
dust attenuation of E(B — V ) derived from SED 
fitting. The top panel of Figure [TO] shows the 
E(B — V) distribution of the highest sSFR galax¬ 
ies (—8.5 < logsSFR/yr _1 < —8.0) and the main 
sequence galaxies (—9.5 < logsSFR/yr _1 < —8.5) 
with solid and dotted histograms, respectively. In 
addition, galaxies are also split into two stellar 
mass bins of 10 10 < _M*/M 0 < 10 10,5 (blue) and 
M* > 1O 1O 5 M 0 (red). The inset shows the cumu¬ 
lative distribution of each sub-sample. The me¬ 
dian values of E(B — V) for the highest sSFR 
bin is 0.45 and 0.55 for 10 10 < M*/M 0 < 10 10 ' 5 
and M* > 1O 1O 5 M 0 galaxies, respectively. How¬ 
ever, those for main sequence galaxies in the 
same stellar mass bins are 0.375 and 0.45, respec- 
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tively. It is clear that high sSFR galaxies are more 
heavily obscured at a given stellar mass and in 
strongly biased regions than th e main sequence 
galaxies. IWhitaker et al. ( 2012h categorized high 
sSFR galaxies as dusty, blue star-forming galaxies, 
and noted that they are possibly merger driven 
starburst galaxies. We also show in Figure [16] the 
distribution of g — J color (bottom) corresponding 
to the rest-frame UV-optical color. The line styles 
are the same as the top panel. This bluer color is 
also confirmed by the best fit SED in Figure [4] 
The bottom panels in Figure |4] are examples of 
high sSFR galaxies. Compared to main sequence 
galaxies (middle panels in Figure[4|, it is clear that 
high sSFR galaxies show significantly higher fluxes 
in the rest-frame UV regime. Although there is a 
significant overlap, high sSFR galaxies are either 
relatively blue g — J color or have a higher extinc¬ 
tion suggesting the presence of an obscured young 
stellar population. 

Although these high sSFR galaxies are most 
likely to be obscured galaxies with efficient star 
formation, their colors could also be strongly af¬ 
fected by the presence of an AGN. From the pre¬ 
vious work based on the clustering, AGN reside in 
dark matter haloes of > 10 13 h~ 1 M Q (e.g., Hickox 
et al. 2009, 2011; Ross et al. 2009b; also see a 
review in Cappelluti et al. 2012 and references 
therein). Interestingly, this halo mass is consis¬ 
tent to our result for high sSFR galaxies. There¬ 
fore, in order to test the possibility of AGN con¬ 
tamination, we check the fraction of AGN in this 
high sSFR bin. Using Spitzer sources mentioned 
above, we select 75 potent ial A GN with IRAC 
colors suggested by IStern et al. ( 20051 ). Among 
these AGN, 15 sources are in the highest sSFR 
bin, and this corresponds to ~8.8% of all IRAC 
sources matched with our galaxies in the same 
sSFR bin. In addition to Spitzer data, we also se¬ 
lect potential AGN in the whole studied area from 
the Wide-field Infrared Survey Expl orer (WISE) 


with a simple color cut suggested bv IStern et al. 


( 20121) . Although the depth of WISE data is not 
enough to detect faint AGN, the wide area is help¬ 
ful to measure the clustering after excluding bright 
AGN. We remove potential WISE AGN from our 
high sSFR galaxies, and then measure the clus¬ 
tering again. The fractions of galaxies excluded 
are 3.6% in the LM bin and 8.9% in the HM bin. 
The re-measured clusterings are identical to those 



E(B-V) 



g-J 

Fig. 16.— E(B—V ) (top) and g—J color (bottom) 
distributions of galaxies with 10 10 < M*/Mq < 
10 10 -5 (blue) an d _M* > 10 10 ' 5 Mq (red). Solid and 
dotted histograms are for the highest sSFR galax¬ 
ies (—8.5 < logsSFR/yr _1 < —8) and main se¬ 
quence ones (—9.5 < logsSFR/yr _1 < —8.5), re¬ 
spectively. The inset shows the cumulative distri¬ 
bution with the same line style to the main panel. 
The highest sSFR galaxies are more heavily ob¬ 
scured and show relatively bluer colors than main 
sequence galaxies at a given stellar mass. 


for the previous measurements in both bins. Since 
the triggering mechanism for AGN is not well un¬ 
derstood yet, we cannot rule out the connection 
between high sSFR galaxies and AGN. However, 
our result for the highest sSFR bin is not biased 
by AGN. 

Above, we noted that galaxies with the highest 
sSFR reside in more massive haloes than main 
sequence galaxies. There are numerous reports 
in the literature that the star formation activ¬ 
ity is enhanced by merger or the tidal inter¬ 


action of gal axies (fM ihos et al. 

2000: Lambas et al. 20031: Ellison et al. 


1992; 


Barton 

2008; 


Hwang et al. EoilJldeue et al. 112012: Kamnczvk et al 


20131: Patton et ahT2013). Daddi et al. (2010t) 
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and lGenzel et al.~ ( 2010h also reported that star- 
burst galaxies driven by mergers form stars more 
efficiently than normal star-forming galaxies. Fur¬ 
thermore, it has been thought that these pro¬ 
cesses more easily occur in dense environments 
such as galaxy groups (McIntosh et al. 2008; 


Perez et al. 


il axy groups (McIntosh ct al. I^UUet 
2009t iPipino et al~ 12014 1. So when 


combined with the observation that these galax¬ 
ies are dustier and form stars efficiently, it sug¬ 
gests that high sSFR galaxies can be generated by 
mergers or interaction. Also, since massive galax¬ 
ies usually reside in massive haloes, high sSFR 
HM galaxies experience the processes in massive 
haloes imprinting the ^0.4 dex high halo mass 
compared to high sSFR LM gala xies as seen in_Fig- 
ure m This is consistent with I Lin et al~l (1201211 
who speculated that galaxies with the highest 
sSFR might be linked to mergers or interactions 
in haloes more massive than those hosting main 
sequence galaxies. Another possible explanation 
for high sSFR galaxies is that they are primarily 
in massive haloes like galaxy with enhance d star 
formation at z ~ 2. iMagliocchetti et al. I ( 2014 1 
pointed out that galaxies forming stars actively 
at z ~ 2 are a different population from similar 
objects at z < 1, based on the comparison of clus¬ 
tering strengths. They reported the halo mass of 
> 1O 13 M 0 , and they evolve into passive galaxies at 
z < 1.5. The halo mass is consistent with that_for 
our hi gh sSFR samples. Moreover, Ponesso et al. I 
( 2015h also find that most IR-luminous galaxies at 
z > 1 reside in the group environment, which is 
also similar to that for our high sSFR galaxies. In 
this sense, our high sSFR galaxies may be in the 
transition phase experiencing a rapid evolution to 
a passive population. 


Now, we recall the HODs in Figure [8] to inter¬ 
pret Figure [15] First, we briefly note the halo 
mass difference measured by the halo model and 
the fit of dark matter clustering, since we directly 
estimate the bias with the dark matter clustering 
for sSFR binned samples. As we discuss in Ap- 
pendixQd] the halo mass from the dark matter clus¬ 
tering fit is between representative halo masses for 
central and satellite galaxies from the halo model. 
If there are more satellites, this becomes closer 
to that for satellites. Furthermore, M e ff of dark 
matter haloes from the halo model is more sig¬ 
nificantly influenced by satellites. In Figure [15] 
the horizontal dashed lines are M e f / values for all 


galaxies in each stellar mass bin. 

For LM galaxies, M e ff for all galaxies from 
the halo model is 1O 12 ' 66 /i _ 1 M 0 which is at the 
higher halo mass regime than the HOD for cen¬ 
tral galaxies in Figure [8] Therefore, if the halo 
mass is larger than this M e ff 1 we are able to 
conclude that satellite galaxies dominantly influ¬ 
ence the clustering. Of LM samples, galaxies with 
— 10 < logsSFR/yr _1 < —9 show halo masses of ~ 
10 13 H~ 1 Mq. This means that these low sSFR LM 
galaxies mainly consist of satellites. In contrast, 
the halo mass for -8.85 < logsSFR/yr _1 < -8.5 
galaxies corresponding to the upper part of the 
main sequence is below M e ff , which implies that 
the significant fraction of them is central galaxies. 

However, in the case of HM galaxies, it is more 
complex than LM galaxies. From Figure 0 it 
can be seen that the central HOD spans a wide 
halo mass range. Moreover, the M e // value for 
all HM galaxies is lO 12 - 83 /i _1 Af 0 , and this is just 
above the peak of the central HOD. Therefore, 
we can reasonably expect a significant contribu¬ 
tion from central galaxies to the clustering. From 
Figure [El the overall trend for HM star-forming 
galaxies is similar to that for LM star-forming 
ones. However, all halo masses do not exceed 
the range of the central HOD of the middle panel 
in Figure 0 and correspond to halo masses of 
N c (M ha i 0 ) > N s (M ha i 0 ), except the highest sSFR 
bin. This also confirms that the important contri¬ 
bution of central galaxies to the clustering. What 
we also find is that the lowest halo mass appears 
in —9.25 < logsSFR/yr _1 < —8.85, which is the 
lower sSFR bin than that for LM galaxies, and cor¬ 
responds to the upper envelope of the star-forming 
main sequence in Figure 0 

Another striking feature is the halo mass of pas¬ 
sive galaxies. We already discussed the depen¬ 
dence on stellar masses in the previous section. 
However, the halo mass of passive galaxies is also 
similar to M e f f of all galaxies in each stellar mass 
bin. This implies that the influence of central pas¬ 
sive galaxies is also significant as well as satellites, 
even in the LM bin. Moreover, the halo masses 
for both stellar mass bins are within the range of 
central HODs in Figure 0 In fact, the fractions of 
passive galaxies are 24% and 50% in the LM and 
HM bins, respectively. Even if we assume that 
all satellites are passive galaxies, these values are 
much larger than the satellite fraction for all galax- 
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ies from the HOD in Figure [9l and confirm the 
substantial contribution of central passive g alax¬ 
i es to the clustering. Additionally, iKrause et, al. 
( 20131) found that red central galaxies tend to be 
in slightly more massive haloes than blue cen¬ 
tral galaxies based on SDSS data. Therefore, the 
higher halo mass for passive galaxies than main 
sequence ones can be explained by the high mass 
of haloes hosting passive central galaxies as well 
as passive satellites. 


The complexity of Figure H5l highlights the po¬ 
tential of HOD analysis to identify subtle envi¬ 
ronmental and evolutionary trends if a sufficiently 
large and well constrained sample of galaxies is 
studied. Based on our results for M* > 1O 1O M 0 
galaxies at z ~ 1, the clustering of main sequence 
galaxies, which shows the lowest halo mass, is most 
significantly affected by central galaxies. Given 
the hierarchical growth of structure and the ex¬ 
pected evolution of main sequence star-forming 
galaxies, it would be expected that these cen¬ 
tral galaxies could become satellite galaxies via 
accretion on to massive haloes or remain as cen¬ 
tral galaxies with the assembly of stellar and halo 
masses. During this accretion, and whilst orbit¬ 
ing as satellites, they would undergo environmen¬ 
tal quenching of their star formation (e.g., from 
ram pressure stripping, tidal harassment or star¬ 
vation), which would explain the high halo masses 
associated with low sSFR galaxies. On the other 
hand, the population with enhanced star forma¬ 
tion activity are in relatively massive haloes of 
group like environments, and a subdominant pop¬ 
ulation at z ~ 1. In our data, the fraction 
of -8.5 < logsSFR/yr” 1 < —8.0 galaxies with 
M* > 10 1o Mq is ~10% of all galaxies, which is 


previous studies (Lotz 

et al. 2008: Bundv et al. 

2009; Conselice et al. 2009; de Ravel et al. 20091 

Lopez-Saniuan et al. 

2009(). Therefore, galaxy 


mergers play a subdominant role for the evolution 
of star-forming galaxies at this epoch, even if we 
assume all high sSFR galax ies experien ce m ergers 
or interaction. Recently iMoustakas et al.' 1 2013 ) 
also concluded that galaxy mergers are not a dom¬ 
inant source of stellar mass growth at z < 1. Ad¬ 
ditionally, the significant fraction of passive galax¬ 
ies are central ones which weakens the clustering 
strength compared to a satellite dominated sam¬ 
ple. These passive central galaxies may originate 


in star-forming central galaxies passively evolving, 
or galaxy mergers with the central galaxy acceler¬ 
ating the consumption of gas. 


Mustek et al. ( 20131 ) concluded that the secu¬ 


lar processes are the dominant mechanism for the 
evolution of galaxies, which means that galaxies 
evolve from a low M*-high sSFR to a high M*- low 
sSFR regime through star formation within the 
galaxy, and that gal axy mergers p l ay a s ubdomi- 
nant role. Moreover. iTinker et al. ( 2013 ) claimed 
that the stellar mass function of passive central 
galaxies has significantly increased since 2 = 1, es¬ 
pecially at M* < 1O 11 M 0 . Our results show that 
main sequence galaxies are mainly central galax¬ 
ies, and a significant fraction of them evolve into 
low sSFR central galaxies together with the de¬ 
creased star formation activity and the increased 
stellar and halo masses. Therefore, our result also 
supports the suggestion that the bulk of M* > 
10 10 Mq star-forming galaxies at z ~ 1 follow sec¬ 
ular evolution (orange arrow in Figure 1151) sup¬ 
plemented by minor mergers with galaxies fainter 
than our stellar mass limit accounting for the 
bulk of their growth . Th is was also suggest ed by 
Noeske et al. ( 2007 ) and Peng et al. ( 20 id) . Ad¬ 
ditionally, asimilai_trend_wasJound at z ~ 2 
( Lin et al. 20121: Sato et aL~| 2014 ]_ and the lo¬ 


cal universe (|Li et al. 120081 : iHeinis et al. I l2009l) 

which implies that this is the main evolutionary 
mechanism of star-forming galaxies over the last 
10 Gyr. Importantly, the magnitude limit of our 
survey allows us to detect both galaxies that would 
give rise to a major merger within any halo so in 
principle it would be possible to constrain the ma¬ 
jor merger rate once the redshifts of a representa¬ 
tive sub-sample of pairs is determined. 


7. SUMMARY AND CONCLUSION 

In this work, we have used deep and wide 
datasets based on UKIDSS DXS and CFHTLS- 
Wide surveys to investigate the dependence of 
galaxy clustering on intrinsic properties and how 
galaxies are linked with their host haloes. The 
main results are summarized as follows; 

1. Using deep and wide near-IR/optical imag¬ 
ing data of the SA22 field, we have con¬ 
structed a mass-limited sample of galaxies 
at 0.8 < z < 1.2. The redshift selection is 
based on photometric redshifts, and galaxy 
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properties such as stellar masses and sSFRs 
were derived using SED fitting. In total, 
this sample consists of 66,864 galaxies with 
M* > 10 10 Mq and logsSFR/yr _1 < —8 in 
this redshift range. 

2. Splitting the selected galaxies at different 
stellar mass thresholds, we measured the an¬ 
gular two-point correlation function and per¬ 
formed the halo modeling to link galaxies 
with host dark matter haloes. We found that 
more massive galaxies reside in more massive 
haloes, and tend to be central galaxies. 

3. The HODs for stellar mass binned galax¬ 
ies were calculated by those for stellar mass 
threshold samples. In all bins, satellites 
are predominantly in > 10 13 /i _1 Mq haloes. 
Moreover, the mass of haloes hosting central 
galaxies is higher for massive galaxies with a 
broader distribution in halo mass than that 
found for less massive galaxies. 

4. The HODs for stellar mass binned sam¬ 
ples were used to calculate the stellar mass 
to halo mass ratio for central and satellite 
galaxies separately. For central galaxies this 
ratio shows a peak at ~ 10 12 ft. _1 Mg that 
drops sharply above and below this halo 
mass, indicating the most efficient stellar 
mass growth at this peak halo mass. On the 
other hand, satellite galaxies significantly 
contribute to the total stellar mass in group 
and cluster environments. 

5. We find an anti-correlation between bias 
and sSFR for low sSFR star-forming galax¬ 
ies that are at or below the main sequence 
(—10 < logsSFR/yr -1 < -8.5), implying 
that star-forming galaxies at around the 
main sequence tend to live in a less mas¬ 
sive halo (Mhaio ~ 10 12 ' 5 /i _1 M 0 ) while 
low sSFR or passive galaxies are more 
likely to be in massive haloes (Afhaio > 
10 12 ' 5 Ii~ 1 Mq). However, we also see a rever¬ 
sal of this relation for galaxies in the highest 
sSFR bin (-8.5 < logsSFR/yi-- 1 < -8) 
that galaxies with the highest sSFRs are in 
dense environments. This can be seen re¬ 
gardless of the stellar mass of galaxies. 

6. We speculate that the bulk of galaxies 
evolves from on or above the main sequence 


of star-forming galaxies to a lower sSFR 
regime as their mass assembles through 
forming new stars and minor mergers. Addi¬ 
tionally, major mergers happen in relatively 
massive haloes, and contribute to the galaxy 
evolution sub-dominantly at z ~ 1. 

Our results are derived from deep and wide mul- 
tiwavelengtlr datasets. Nevertheless, this work is 
based on the only photometric data for relatively 
massive galaxies in a specific redshift bin and thus 
it is difficult to avoid contamination. In the future, 
deeper and wider datasets such as those from the 
UKIDSS, VISTA, Subaru Hyper-Suprime Cam¬ 
era, Pan-STARRS and LSST surveys will provide 
an opportunity to investigate more details about 
the relationship between various galaxy properties 
and their host dark matter haloes with less mas¬ 
sive galaxies and in various redshift bins. Addi¬ 
tionally, spectroscopic surveys will also allow fur¬ 
ther progress on the clustering analysis. In terms 
of parameterized HODs, we have assumed that 
the central HOD becomes unity after a certain 
halo mass. However, AGN feedback may change 
the shape of central HODs, especially the maxi¬ 
mum mean number for central galaxies. Therefore 
it may also be necessary to modify the standard 
HOD work in order to directly compare observa¬ 
tions and models once the samples have increased 
sufficiently. 
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Table 2: Results for sSFR binned galaxies with 
stellar mass thresholds. Column (1) is the stellar 
mass threshold in Mg, column (2) shows the sSFR 
range with a logarithmic scale in yr -1 and column 
(3) is the number of galaxies in each bin. Column 
(4-5) are the measured bias and the estimated halo 
mass, respectively. The unit of halo masses is in 
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Table 3: The same table with Table 0 but for stel¬ 
lar mass binned galaxies. Column (1) represents 
the stellar mass range. 


M, 

sSFR range 

Ngal 

bias 

i^halo 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

10 lu < M, < lO 1 " 5 

-8.50 - -8.00 

5082 

i-9o7o:o9 

12.6941°;°?; 

(LM) 

-8.85 - -8.50 

6675 

1 yq+o.08 
x - ‘ 0.09 

12.5851°;°® 


-9.25 - -8.85 

6100 

2.56l°;; l i l 

13.1491°;°® 


-9.70 - -9.25 

4049 

2.50lS;“ 

13.1141°;°?° 


-10.1 - -9.70 

2394 

9 QQ+0.17 
z,.»O _ 0 . 1 8 

13.3291°;°?° 


<-10.5 

9095 

1 Qc;+0.07 

1.00_Q Qg 

12.64318;°!? 

o 

o 

A 

A 

o 

-8.50 - -8.00 

1257 

9 A A +°- 18 

13.0811S;!;? 

(HM) 

-8.85 - -8.50 

2418 

9 1 Q+0.14 

12 877 + 0.110 
lz -° 1 1 -0.103 


-9.25 - -8.85 

2807 

i-94±g;l! 

1 9 7OCC+0.119 
1 Z. 1 — 0.110 


-9.70 - -9.25 

2260 

2.04l°;“ 

12 . 80618 41? 


-10.1 - -9.70 

1527 

9 0^+919 
Z.ZO- 0.21 

12.9617841? 


<-10.5 

12280 

2 ia+0.06 

12.85T78;8J1 


This 2-column preprint was prepared with the AAS lAT^X 
macros v5.2. 
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Table 1: The HOD parameters for stellar mass threshold galaxies at 0.8 < z < 1.2. Column (1) represents 
the stellar mass threshold for each sub-sample in Mq, columns (2-5) are the best fit HOD parameters and 
column (6) is the number density of galaxies in 10 -4 /i 3 Mpc -3 . Columns (7-9) show the derived quantities 
base on equations (5)-(7). The final column is the quality of the HOD fit in terms of % 2 per degree of 
freedom. All dark matter halo masses are in h~ 1 M Q with a logarithmic scale. 
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A. HOD MODELING WITHOUT n g CONSTRAINT 


As desribed in § 13.31 the halo model with three free parameters was applied to derive the property 
of haloes hosting our galaxy sample. In this case, M cut was determined by matching the observed galaxy 
number density with given parameters. In order to check the influence of this constraint to the best fit result, 
we investigate the best fit parameters without this constraint. Therefore, we perform the halo modeling for 
M* > 10 lo M e and M* > 1O 11 M 0 galaxies with four free parameters ( a cut , M cut , M 0 and a). In this case, 
rig is derived by the best fit HOD parameters describing the clustering only. 

Figure [T7] shows the angular correlation function (upper) and the HOD (lower) for each galaxy sample. 
All symbols are the same as Figure [5] However, the solid line in the upper panel and the HOD in the lower 
panel are the result based on the halo model with four free parameters. The best fit parameters are listed 
in Table [3J Comparing the values in Table [T] and SI the derived values of b g and M e f f are identical, and 
other halo mass parameters also show the same trend. In addition, the central HOD for massive galaxies 
(a cut =0.4) still shows the gentler shape than that for low mass galaxies ( a cut =0.1), which indicates that 
the constraint by the galaxy number density does not affect the trend of the fitting results. 


However, M cut and Mo without the n g constraint are approximately a factor of 2 smaller than those 
with the constraint, and <j cut values also decrease from 0.5 and 0.6 to 0.1 and 0.4 for M* > 1O 1O M 0 
and M* > 10 n M 0 , respectively. In addition, the number densities calculated by t he m odel without the 
cons traint are higher t han observed ones. This mismatch was already reported in Matsuoka et al. ( 201 1;) 


and Wake et al. ( 201lh . Unfortunately, it is not yet answered what is the main reason leading to this 


mismatch. More studies are necessary to resolve this problem. 




1 M,> 10 10 0 M 0 ! 

M*> 10 110 M 0 e 

- r\s 






0.001 0.01 0.1 10 . 0010.01 0.1 1 


0(deg) 



M halo (h- 1 M 0 ) 


Fig. 17.- The same plot with Figure [6] for M * > 1O 1O M 0 and M* > lO n Af 0 galaxies. All symbols are the 
same as Figure [6] However, we used four free parameters for the halo modeling. As a result, the trend of 
derived parameters depending on M* is same to Figure[6l although values are different. 
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B. HALO MODELING vs. FITTING DARK MATTER CLUSTERING 


In order to link galaxies to their host dark matter haloes, we performed the halo modeling for stellar mass 
threshold samples and the direct fit of correlation functions of dark matters for stellar mass binned samples 
and sSFR binned samples. Here, we apply the later method to the stellar mass threshold samples, and then 
compare bias and halo mass from this fit (quoted as DM fit, hereafter) to the best fit result from the halo 
model. This comparison will provide a guideline for our analysis. 

Figure [18] shows the comparison between the halo model (open symbols) in Table [I] and the fit of dark 
matter clustering (filled symbols). First, we are able to notice that the bias is consistent independently of the 
method (upper panel). However, in the case of halo masses, they show different results. The halo mass from 
the DM fit is close to M e ff from the halo model, but shows a discrepancy, especially for low mass galaxies. 
This is easily explained by the fraction of satellites, since satellites are in massive haloes. In addition, M e ff 
is close to Mq at the low stellar mass regime, but M cut at the high mass regime, which is also explained by 
the same reason. 

Although the bias estimated by fitting the dark matter clustering do not perfectly represent the mass of 
haloes hosting galaxies selected, representative halo masses for central and satellite galaxies well bracket the 
halo mass by the DM fit. Therefore, the halo mass for central galaxies is always lower than that by the DM 
fit. 



10 10.5 11 

log(M./M 0 ) 


Fig. 18.— Comparison of results from the halo model (open symbols) and the fit of dark matter correlation 
function (filled symbols) for stellar mass threshold samples. Upper and lower panels are for bias and halo 
mass, respectively. For the display purpose, the errors from the halo model are not displayed (see Table |T] 
for these values). 
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Table 4: The HOD parameters for stellar mass threshold galaxies. However, the halo modeling is performed 
with four free parameters without the constraint by the galaxy number density. All units are the same to 
table m 
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